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vaux modernes, relatifs aux mêmes questions, il faut mentionner 
le Mémoire queEncke a rédigé en i83o, d'après un cours de Gauss 
auquel il avait assisté en 1812, et qui vient d'être réimprimé, 
ainsi que le Mémoire posthume de Gauss sur l'interpolation; 
puis, surtout, les recherches de Cauchj (1 835- 1 853) et celles de 
M. Tchébychef (1854-1875). Le Rapport de M. Merrifield, pré- 
senté à l'Association britannique en 1880 (*), n'en fait point 
mention. 

I. 

1. Considérons les deux séries parallèles 

3^0 *^l *^2 • • • 

7o 7i 72 ••• 

et désignons par Dj^o? D^< ? • • • , D^j/q, D^y^ , . . . , ou plus simple- 
ment par Do, Di, . . . , Djj, D^, ... les rapports 

n _ yi — .n n _ yt—n 

^1 — Xo Xi — Xi 

/x^ ;n2-^lZL^ n2 D2-D1 

(i) { ^0= ' ^1 = > • • • > 

X2 — Xq X-^ Xi 

373 — Xq Xi, — Xi 

dont l'expression générale 






(2) 

a pour dénominateur la différence des valeurs extrêmes de x qui 
ont encore concouru à la formation de chaque rapport, comme le 
montre le Tableau 

^0 yo 

Do 

^1 JKi D? 

D, D3 

Xi 72 Df DJ 

D2 DJ 

I>3 

(') Report on tlie présent state 0/ knowledge 0/ the application 0/ qua- 
dratures and interpolation to actual data. 
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En l'écrivant sous cette forme 



• > 



il saute aux yeux que noire algorithme ne fait que traduire, ou 
résoudre en ses éléments, la formule d'interpolation 

qui devient 

\ y = yo-h(x — xo)Do-h(x — xo)ix — xx)D^ 
\ -h{X — Xo){x — Xi){x — X2)T)l-\-.,,, 

si nous mettons ^, y, ^o» J^o? - • • à la place de Xq, y^, x^^ y^^ . . . , 
en reculant les indices d'une unité. C'est la formule générale d'in- 
terpolation de Newton, à peu près telle qu'on peut la tirer du 
5® lemme du Livre III des Principes {Invenire lineam curvam 
generis parabolici, quœ per data quotcunque puncta transit,) 
On peut aussi commencer par poser 



d'où 



yt= DS-+-(^t — ^o)i>J-H(^2 — ^o)(a?2 — a7i)D5, 

de sorte que la signification des constantes DJ résulte du rôle 
même qu'elles jouent dans la formule. Mais nous allons les consi- 
dérer encore sous un autre aspect. 

:2. On voit sans peine, en développant, l'une après l'autre, les 
expressions (i) des quantités DJJ, qu'on aura 

Uo = ' > 

^0 — ^1 ^1 — ^0 

D» = ^ : + . ^ + - ^' 



et, en général, 

(4) I>ï=7 — ^'"-^ —t ^^■' 

(Xo — X\) , . . (J^o — ^n) t««^/t — Xq) . . . {Xfi — X,i—i) 

Le rapport D" est donc une fonction symétrique des quantités 
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Xqj . . . , Xnt yo-i • • • > yn-i dont nous sommes libres de permuter 
les indices. C'est ce qui justifie la nouvelle notation 

Do = (oi), bi = (i2), ..., DJ=(oi2), Dî = (i23), 

analogue à celle d'Encke dans son Mémoire de i83o, et que fera 
mieux comprendre le petit Tableau ci-après : 

o 

(01) 

I (oi'i) 

(12) (0128) 

a (128) (01234) 

(28) (1234) 

3 (234 (12345) 

En écrivant o, i , . • • , pour Xq^ x^^ . . • , on a, par exemple, 

^. (012) — (128) (012) — (oi3) (028) — (128) 

(0128)= 5 = 5 = = 

' o — 3 2 — 8 o — I 

Si l'on suppose maintenant que y est une fonction entière de x 

y =Ao-+-Aia7 -+- A2 rr^ -f- . . . , 

il est facile de voir que DJ représente la partie entière du déve- 
loppement de l'un quelconque des termes de l'expression (4); on 
aura, par exemple, 

(5) DJ=E '^' 



{Xq — Xx) . . . {Xq — Xn ) 

OU bien, en développant, 

DJ = E(A„-h A„4.ia7o-H...)(n- — H-...)...(n- — -4-...), 
DJ = Art -H. . .H- Ap(a7o, x^t . . ., ar^ )''-" + . . .. 

En appliquant l'opération D'' aux termes successifs de j^o? on 
voit, en effet, que A^ a en facteur le poljnôme 

(6) D«a?J = (^o, ar,, .,.,Xn)P-" 

qui se déduit de {xq + ^^ +. . . + Xn)P~" en y remplaçant tous 



— s — . 

les coefficients par l'unité. On a, en particulier, 

D«ar;+* = a:o-^.^l-4-...-^ir„, 

Les D s'annulent si leur exposant dépasse celui de x. Si doncj^ 
n'est que du degré p, nous aurons DP'^^yo = o. 
Supposons que y soit du degré /i, 

(7) y = ko-h A^ix-h,. ,-\- AnX^' 

En ajoutant à cette relation celles qui s'en déduisent en remplaçant 
Xy y par Xq^ y^^ . . . , Xni yni nous avons n 4- 2 équations qui per- 
mettent d'éliminer les n + i coefficients Ao, . . . , A;,, et le résultat 
de l'élimination, c'est l'équation 

(8) D«+i^=o 

qui, développée à l'aide de la formule (4), devient 

y 

{X — Xq){X — Xx),.,{X — *Xn) 

^ To +...-+- Zl =0, 

{Xçs — X){x^i — Xx),, .{XQ — Xa) '" {3rn~X){Xn — Xo).. .{Xn — Xa-i) 

ou bien 

(X — Xi) ,. .(X — Xfi) 



(9) y=ro 



{Xq Xi) . . . (Xq — X,i) 



C'est la formule d'interpolation de Lagrange, qui exprime y en 
fonction de x et des quantités données Xq, • • • , Xn, yo, • • • ? J^/i- 
Lagrange Ta donnée en 1 790, dans ses Leçons élémentaires (Jour- 
nal de VEc, PoL, t. II, p. 277). D'après Boole, on la trouve aussi 
dans Euler (Z?^ eximio usu methodi interpolationum in sérier um 
doctrina). 

Nous avons, d'autre part, 

Yiny _ ^^ -_ const., 
et, par suite, 

(10) D«7 = D«7o, 
équation qui, développée, devient 

(11) |^=^o-^(^ — ^o)^^o 
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C'est Ja formule d'interpolation de Newton, que nous avons déjà 
démontrée plus haut [formule (3)]. 

3. Il faut mentionner ici une Note de M. Emorj M*^ Clintock, 
publiée dans le tome II de Y American Journal of Mathema- 
tics (1879), sous ce titre : A new gênerai method of interpola- 
tion. L'innovation consiste simplement à introduire, au lieu des 
rapports (1), les suivants : 

Xx — ocq "^ x^ — a7o 
D* = > 1^3 = > •••> 

^ X^ — Xx Xz— Xi 

dont l'expression générale est 

(12) Dr' = ^^^ — —■ 

Xi — Xn 

On a ensuite 

y = yo-\-{x — xq)D\-{-{x — xçi){x — xx)iyi H- 

Ce mode de calcul rend, selon l'auteur, la besogne plus mécanique 
(work more mec/ianical). On constate aisément que les têtes de 
séries D], D^, D^, ... de M. M*^ Clintock peuvent être repré- 
sentées par le second membre de la formule (4), et coïncident, 
en valeur, avec les rapports que nous avons désignés plus haut 
par Do, D^, DJ, 

4. On sait, d'autre part, que la formule de Lagrange a été géné- 
ralisée par Cauchy. En cherchant à exprimer y^ non plus par une 
fonction entière,' mais par une fonction rationnelle de x, de la 
forme 



03) 

il trouve 



a -h bx -h. . .-h /ix^^ 

a-f- ^x -h. .,-\- sa:"*' 

yo.Vl • • • YniK^ — ^m 4-1) • • • (^ — X m->rn ) 



__ ('^0 X m-\-\ ) • • • ( ^0 ^m + n ) • • • i^fn ^/«-i-t ) • • • ( ^m — ^m-^n ) 

""* 7ori " -yin-xio^Q — X), ,.{x,n-i — x) 

\Xq ^m) • • • (^0 •^/«-+-/t)« • • K^m—i Xf/i) • • • {^m—i — ^m+n) 

Cette formule a fait l'objet d'un Mémoire de Brassinne(/o;//72<7/ 
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de Lioiiville, t. XI, 1846); elle a été aussi étudiée par Jacobi 
{Journal de Crelle, t. XXX, i845) et par Roscnliain {ibid.). 

Je me bornerai ici aux remarques suivantes. Si l'on veut utiliser 
l'expression (i3) pour l'interpolation, il convient de mettre le 
numérateur et le dénominateur sous la forme adoptée par Newton. 
Posons 

Yi =^1 -h(a7 — iri)Di-+-(ir — Xi){x — Xi)D\ -^ . . . , 

•• • > 

et 

a(r — Yo)-f- p(^— ^o)(r — Yi)-»-ï(^ — ^o)(a7 — ir,)(7 — Yî)4-... = o, 
d'où 



y ^ 

OL-\-{x — X^)^ -\-(x — Xq){X — X^Y(-<r,,. 



OU bien 



_ ajo-t-Çr— yo)(gDo-hP7i)-i-(3r — jro)(3: — .ri)(gDi;-+-pDi-4-Y.r8)-^-" 
•^ a->r{x — Xq)^ -\-{x — Xq{x — d7i)Y-i-... 

En faisant m -\- n =p^ le nombre de valeurs de y dont on aura 
besoin serait? + i. Si alors on fait successivement^ = Xq, ^o . . ., 
Xpj le premier membre de l'équation 

a(7 — Yo)-hP(a7 — aro)(y — Y,)H-...= o 

se réduit à o pour x = Xq^ . . ., Xn, parce que Yq, Y<, . . . , Y,„ 
(qui sont des degrés /i, n — i, . . . , n — m) coïncident avec les 
valeurs correspondantes dey ; mais les substitutions x = x^^^ ? • • • > 
Xp donnent les conditions 

aDJ+» -4- pDï H- vD?-i H-. . . = 0, 

aDJ -f- pDf-» H- yDÇ-2 -h. . . = o, 

et Ton a ainsi m + i équations entre lesquelles on peut éliminer a, 
[3, Y, .... On voit que le dénominateur de y s'obtiendra en for- 
mant le déterminant 



I 


(X — Xi^) (X ■—Xo)(x — Xi) 


i>{; 


J)P 1 \)n-2 


l>'o ' 


\)P 2 i),>-:i 



- 8 — 
et le numérateur en y remplaçant la première ligne par celle-ci 

La formule ainsi construite pourra être, en quelque sorte, ren- 



I I 



versée, en remplaçant partout j',)©, . . . par -, — ? • • •; transfor- 
mation qui peut aider à la simplifier. On trouve ainsi, par exemple, 
pour m^= n = i, avec trois valeurs données de^. 



_ roDi-4-(:c — aro)(D|Do— ytP;) 

J — —— 

OU bien 



y D,-(^-;ro)DJ 






S. Dans le cas où ^0» ^\^ • • • forment une série arithmétique, 
si nous prenons l'intervalle constant (;r, — ^0) pour unité, nous 
aurons, en désignant par A, A^, ... les difTérences successives, 



I . '2 . . . n 
En faisant x — Xo= t, la formule de Newton devient 

(i4) jr, = 7o+<Aro+ \,.^^ "7o-f-- ;-;^^;3 ^ A3jo 

OU bien 

où t^ ^2> ^3? •• • sont les coefficients du binôme. On suppose ici 
^>^oj ^>o. Pour :r <; ^To, on fait ^0 — x=^t^ et la formule 
devient (en remplaçant t par — /) 

(I5) y.,=y,-tt^y^-\- \ t^^y,--^ -^ ^A3jo4-..., 

ou bien 

7- r = ^0 — ^ A H- ( f -f- 1 ) î A 5 — ( ^ 4- 2 ) 3 A J -+- . . . , 

OU encore, en avançant les indices, 

jo=y^— ^ Af-+-(<-f-i)2A| — (f -h2)3A?-H 

On peut donner une remarquable transparence à toutes les for- 
mules de cette théorie par l'emploi de la notation symbolique 
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de Leibniz, que Lagrange a développée dans son Mémoire de 1772 
{CEuvreSy t. III, p. 440> ^^ 4"^ X^^vià à devenir d'un usage cou- 
rant {voir BooLE, Calculas of Jinite différences). Dans cette 
notation, 

j, = (H-A)7o, 

et la formule de Newton s'écrit 

(16) 7^ = (i-+-Ayjo, 

ou bien, si nous omettons y, en désignant j^o par ^X? comme le fait 
Lagrange, 

On aurait de même 

(17) r-/ = (i + -^)-'ro, jKo = (i + A)-0'^. 

Les différences Aq, Ajj, . . . , dans ces formules, sont les différences 
successives qui, dans le Tableau, descendent en diagonale, à 
partir de jKo^ 

y\ •" ^0 

72 ... Af 
... A2 

« 

Dans les formules (i4) et (i5), elles servent respectivement à 
interpoler en avant et en arrière. Si, pour interpoler en arrière, 
on voulait suivre une marche rétrograde, symétrique de la marche 
progressive, il faudrait employer les différences qui montent en 
diagonale (en changeant le signe des différences d'ordre impair), 
et les formules (i4) et (i5) s'écriraient 

'(^ — .• t(t — l)(t — 2) ,, 

(,8) 7-.=7o-M->-H-i^Al,^ ^^^^3 ^ Ai3 
(19) yt = y^-\'t\-i-h \ ^ ^ AJa-f- -^^ — —^ iM, 



1.2 



1.2.3 "' 



ces transformations s'obtiennent, sans calcul, par l'inspection du 
Tableau des différences. 
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6. Reprenons la formule (3), où les intervalles sont quelconques. 
En écrivant a, b^ c^ ... à la place de Xq, x^ , x^^ . . . , et faisant 
usage de la nouvelle notation des rapports D'', qui a été expliquée 
au n" 2, elle devient 



y = ya-^{^ — a){ab)->r{x — a){x — b){abc) 
-\-{x — a){x — b){x — c){abcd) 



(20) 



Il est clair que les lettres a^ b. c, ... pourront représenter les 
valeurs données de x dans un ordre arbitrairement choisi. Si nous 
écrivons a, 6, a', 6', ... au lieu de a, bj c, d, . . ., en faisant 
coïncider la suite naturelle des éléments ^ avec l'un ou l'autre des 
deux groupements symétriques 

... a" a' a b b' b" . . . , 
... b" b' b a a' a" . . . , 

il est visible que la formule, qui devient alors 



{ y = j^-+-(ip_ a){ab)-\-{x — a){x — b){a'ab) 
\ -i-(x — a')(x — a)(x — b){a abb')-h. . .^ 



(21) 



n'introduit (par une marche en zigzag) que les différences placées 
immédiatement au-dessus ou au-dessous de la ligne horizontale, 
menée par le milieu de l'intervalle ab : 



a' 



a 



(a'a) 

... (a ab) 
(ab) ... (a'abb') 

... (abb') 
{bb') 



b' 



b' 



b 



a 



a' 



{bb') 

{abb') 
(ab) ... (a'abb') 

{a' ab) 
{aa') 



Cette formule convient surtout lorsque x tombe entre a et b. 

Au lieu de varier les groupements, on pourrait aussi en con- 
server un seul, le premier, et faire, en revanche, des permutations 
de lettres dans la formule (21). Les deux permutations 



a" a' a b h' b" 
b' h a a a" a'" 
b" b' b a a' a" 
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donneraient respectivement 

(2-2) I y = ya-i-(^ — a){aa')-h(x — a)(^ — a){a' ab) 

( -h(x — a'){x — a)(x — b)(a''a'ab)-^,. ., 

(-23) ^ y=^yb-^{x—b){ab)-\-{x — a){x^b){abb') 

^ ^ \ -h(x — a)(x — b)(x — b')(a'abb')-h.... 

En combinant (21) d'abord avec (22), puis aver (23), on trouve- 
rait successivement 

\y=ya^{x — a) '- ^ '- -^{x—a){x \{a'ab) 

('lA) ^ ^ \ ^- / 
) / ,x/ ,r ,Aa'abb')-^(a''a'ab) 
\ -h{x — a){x—'a){x — b) ^^ -h..., 

( ya-^Vh I a-^b\, ,. , ., , (a'ab^-^(abb') 
\y= \ -^y" —j{ab)M^-a){x^b) 

I -^{x — a){x— b)lx ; j(a'abb')-h 

s 

7. L'aspect de ces formules se simpb'fie beaucoup lorsque les 
valeurs données dey sont équidistantes, c'est-à-dire séparées par 
des intervalles égaux de l'argument x. En faisant x — a=dtzl, 
b — a = i, et confondant a avec Xq, les formules (21) et (22) 
deviennent 

y,=y,^t^,^-—-l,l,^ —.^ A3, 

(.6) <' ^(.-H,)/(^-.)(.-.)^ ^^^ 

1.2.3.4 



et 

Al, H 

1.2 1.2.3 



yt =7o-h « A_iH -' Al, -+- ^^ -^. -' Alj 



(^7) < + rïX4 ^s-^---' 

j._. = yo~^A.,-^_^Al,-^--;-3--^Ai,4-.. 

De même, (3.3) nous donne 

{ ,. X. '(' — '">»% '(' — oc — ■>-) ., 

|r. = r. + («-')Ao+-^^A5 + ^3 Ai, 
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Au lieu d'établir ces formules directement par l'interprétation des 
groupements, on pourrait les déduire des formules (i4) et (i5) 
par substitutions successives» en faisant 

et ainsi de suite ; mais la marche que nous avons suivie, avec Gauss, 
est plus simple.. 

Je ne citerai que pour mémoire la formule de Laplace {Théorie 
analytique des Probabilités, p. i3), qui peut s'écrire 

r^ = 7o-H^A_;.H -— . ^A_22,.H ■ — f^T -'A23,.-H... 



OU bi 



len 

t(t — ar — i) ,, 

yt=ro-^t^r-^ Al^-H.... 

Par la combinaison des formules (26) et ('2^), on obtient la 
suivante, qui correspond à ('24), 

>A a "T" • • • • 



\ I . '2 . 3 '1 1.2.3.4 

C'est la formule dite de Stirling, mais qui était déjà connue de 
Newton (Prop. lll de la Méthode différentielle). Elle est surtout 
propre à l'interpolalion entre un nombre impair de quantités 
équidistantes. Pour un nombre pair, on aurait la suivante, qui 
résulte de la combinaison de (26) et (28), et qui correspond 

à (25) 

j^,= Zi±i:iHH(.-i)Ao 

I ,1 2 1.2.3 

Cette formule, que Ton trouve aussi dans Stirling, peut s'écrire, en 
simplifiant les deux premiers termes, 

ni) r -V ^t\ • fUzii)^l±^i . 
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ou bien encore, en comptant t à partir du milieu de Tinlervalle 

\ I . !i . 1.2.3.4 '2 

On trouve, dans le Recueil de Tables astronomiques de Peters, 
deux Tables qui facilitent l'application des formules (29) et (3i). 
Pour la formule de Newton, lesvaleurs des coefficients sont données 
par Vega, Barlow, etc. Voir aussi Woolhouse, On interpolation 
(Londres, i865). 

En faisant ^ = o, (32) donne la formule connue 

(33) ^.= Z5L±Zl_I ^i±^+ 3 -^-1-+-^-^ 



* 2 82 128 2 

8. L'avantage des formules qui précèdent, c'est qu'elles mettent 
en jeu des différences « symétriques », c'est-à-dire échelonnées le 
long d'une ligne de symétrie (d'une ligne horizontale) du Tableau, 
ce qui favorise la convergence. Mais, lorsqu'on s'arrête aux for- 
mules (26), (27), la symétrie n'est pas apparente dans les notations, 
et la marche est alternante : pour passer de Aq à A^,, l'opération 

est -> et de A^, à A^,, simplement A. Sur la même ligne, rigou- 
reusement, on ne peut passer que de Aq à A^, , par l'opération , 

qui ne se dédouble pas : 

A . A3 

• • • «lA I • • • *Jk A • • • 

^0 • • • ^-1 • • • ^-S 

A« A3 

Pour retrouver la symétrie, on a recours aux différences que 
nous appellerons médianes (c'est-à-dire placées sur la ligne mé- 
diane du Tableau), en comblant les lacunes par les moyennes 
arithmétiques des différences de même ordre, placées immédiate- 
ment au-dessus et au-dessous de la ligne, comme le montrent les 



— a — 

formules (29) et (3o). La ligne menée parle travers de y© contient 
les did'érences médianes 

(J4) > -^-1» » -^-j» .••> 

taudis que celle qui passe entre yo ely^ renferme 

' '2 2 2 

Afin d'abréger l'écriture, et de mieux marquer la symétrie de 
l'algorithme de l'interpolation, on introduit souvent, dans la nota- 
tion des différences, l'argument moyen, qui correspond à la ligne 
sur laquelle la différence est placée, au lieu de l'argument de 
départ, qui figure dans l'ancienne notation. Les différences mé- 
dianes, placées sur une même ligne, ont alors toutes le même argu- 
ment. Si, pour plus de clarté, nous les désignons par la lettre D, 
celles qui forment la série (34) s'écriront 

D(o), D2(o), D3(o), D*(o), ..., 

et celles de la série (35) 

Do(i), D(i), D»(|), D3(i), DHi), ... 

ou, si l'on préfère la notation d'Encke, 

/'(o), /'(o), /'"(o), /-(o), ... 
et 

/^(i), /'(i), /^(i), ri\). r\ï)^ .... 

Hansen, dans un Mémoire de i865 {Retationen zwischen 
Summen iind Diffère n zen) met Targument moyen en indice; il 
écrit, par exemple, Aj^«. Cet usage est critiqué par G. Boole, qui 
le déclare anti-scientifique {unscientific and wliolly in défiance 
of convention)', mais la notation est commode et ne causera pas 
de méprises, si l'on a toujours soin de mettre l'argument moyen 
entre deux parenthèses; avec cette notation, les formules (26) 
et (27) s'écrivent 

I r..=^.±<D(±i)+ £^ D.(o)± lilillliiilli} Da(±î)+..., 
(36) < '•'■'•^ 



ou bien 

/(±<)=/(o)±</'(±i)+<, /'(o)=t(< + i),/'"(±;)+..., 

/(±0 = /(o)±</'(Ti) + (« + i)î/'(o)d:(< + i ),/""(+!) 



•» 



La formule de Stirllng s'écrit (pour ^<o) 

(37J.yi=r. + <D(o)+-i^D'(o)+ 'Y,^ 3^ D»(o)+ Y..,.3.4 P'(») 
el la formule (3i) 

Enfin (33) devient 

(39) 7i=Doa)-iD2(A)+-l_Dv(';__l-D6(i) + .... 

La formule de Stirling offre l'avantage de se prêter, avec la 
même facilité, à l'interpolation en avant et en arrière, puisqu'il 
suffît de changer le signe des termes d'ordre pair pour l'interpola- 
tion en arrière. 

En la développant suivant les puissances de t, on trouverait 

(r/=7o-+-«(D-iD3-i-3LD5-...) 
(40) ■+-^î(iD«-^Dv + ï|;jD6_...)4-«3(lD3_J^D5H-...) 

( -H «H ïLD^__j_D6^_...)_^_^5(^D5 — ...)■+•.... 

9. Pour les applications que nous avons en vue, il importe de 
nous rendre compte des relations qui existent entre les différences 
ordinaires A et les différences médianes D, dont le caractère mixte 
s'oppose à l'emploi de l'algorithme habituel. 

D'après les définitions du numéro précédent, nous avons 

Do(o)=Ai>, D(o)= ^^-^^"S 

A», _|_ A3 

D>(o) = AÎ,, D3(o)=-^i— ^i, 



D2«(o)= AIJ, D2«-^i(o) = 



î 



•1 
On aurait aussi, comme il est facile de voir 

j D2«(o)=Al«, D2«+i(o) -l-|D2«H-2(o) = A1S+», 
( D2«(n)=A5", D2«M(^)4_^-D2«-^-2('n) = a;«-^\ 
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et par la formule (17) 

(17) Ai„=A*(H-A)ô^ 

Pour simplifier récrilure, nous supprimerons l'indice o el Targu- 
ment (o) toutes les fois que cela pourra se faire sans nuire à la 
clarté. On trouve alors, en se servant de la formule (17), 

5tD = AH =i-hA — (1-4- A)-», 

1-+- A 

(4a) \ D2= -^ = i-f-A-+-(i-4- A)-> — 9., 

I -h A 

(i-f.A)=^»= iihDH-jD2, A = D-T-tD2. 

Pour une différence quelconque, d'ordre pair ou d'ordre impair, 
on trouve de même 

( D2« = A2«(i -h A)-« = [(i-f- A)2 — (i -h ^)~'^T\ 

(43) < 2Dî'»+i=A2«+i[(i -H A )-'*-+-( 1-4- A )-«-'] 

( =[(i-+-A)-(i-i-A)-i][(i-T-A)^-(i-4-A)"5]*". 

OnadoncD2« = (D-*)", et 

mais non pas (D)^" = D-" ; on ne peut donc associer directement 
deux exposants impairs, ni élever D à une puissance, comme il est 
permis de le faire pour D-, On a ici 

(44) < (D)2'^=(D2^lD'0''= D2/Mi-f-in2)", 

( D.D2«+î= D2«+2-h{D2"-^^ 

En développant les expressions {/\'i), on trouve ensuite 

+ "("-^" A«+î _ n(/t + i)(ft + a) ^,„^3 ^_ 

(45) < 
^'* ^ ^ D2«-M= A*«-»-i — (n-f- J)A2«-«-2 

^^ + '^^ A2/^-^3 _ (n-^i)(n^',)(n-^l) 

., •" i . • • ^ 
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par conséquent 

D =A — lA2-+-|A3— jA^-t-..., 

D2 = Aî— A3-4- A^ — .... 

D3=: A3 — fA*-+-aA5— ..., 

D* = A* — 2A«-h3A6— 

Pour arriver à exprimer les A par les D, nous appliquerons l'opé- 
ralion D^" à la formule de Slirling, en y faisant ^ = /z, ce qui nous 
donne 

D2n(/l)=D2/ï(0)-f-/lD«'»+I(0)H- — D2«+2(o)-4-.... 

Le premier membre étant remplacé par A2"(o), nous pouvons 
supprimer partout Targument o, et il vient 

^ a 2.3 



ou 



bien 



H-(/l4-l)3D*«+3_| ?^_^ (;i_^_2)4D*«+*-h 

2AI -+- 4 

On a, d'autre part, 

et, en tenant compte des relations (44) dans la rencontre des expo- 
sants impairs, on trouve 

+ ^(^ + D2«+3 ^ /t(/l-M)(/l + i) ,),„^, _^^ 
1.2 1.2.3 

OU bien 

A««+«=: D2«+*^- ^^^"^ ' (/l + l)D2«+2 

2 /l -h 2 

-+-( 71-4- 1)2 D^'+s-f- ^^"^ ' (n H-2)3D2«+* + .... 
^ ^ 2 W -H 4 ^ 

Il s^ensuit qu'on aura 

(i -h A)IA = I -h^ D«'»+' ( (Jt2,i-Hi -{- n.iiî„-h n^. fXs/i-î H- . • • ) 



^D2'*-^[([JL — i)2„_,-i-/î(a — l)2„_2-t-...]. 



\{. 
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Or on a généralement (en désignant toujours par rtp les coefficients 
de la puissance /i'''™*' du binôme) 

« 

Il s'ensuit que 

1 

ou bien 



Or ie% formules (42) montrent qu'on aura (1 + A)~* en rempla- 
çant D par — D dans (i + A); on trouve donc 

( 43) <« ^- A)*H' — I = V (±1 Dî^+i 4- ^^^ Dî^+îVjjL -H /i)s«+i. 



Cette relation nous sera utile dans la suite. 

10. La méthode de Mouton que Lagrange développe dans son 
Mémoire de 1792 (*) consiste à déterminer, par les différences de 
la série donnée, celles de la série interpolée, et à remonter de 
celles-ci aux termes cherchés, de sorte que le calcul se réduit à 
des additions successives. 

Les formules s'établissent très facilement à l'aide de la notation 
s}rmbolique. Étant donnée une série de termes équidistants Tq, 
T|, •.., il s'agit d^tnterpoler les termes également équidistants 
t^y t^y . . . , de façon à subdiviser l'intervalle primitif en m inter- 
valles plus petits. Nous aurons 

et, en désignant par A les différences des T, par 8 celles des /, 

par conséquent 

'(49) i-f- A =(i-l-o)'«. 

(') Œuvres, t. V, p. 663. Voir aussi le iMémoire de 1778 5m/' les interpolations 
{Œuvres, t. VII, p. 535) et le Mémoire de 1783 Sur une méthode particulière 
d^ approximation (t. V, p. 517). 
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Il en résulte que 

et aussi 

A* = [(i -f- 8 )'« — i]* = (m 8 4- m, 8« -+-. . .y. 

On trouve ainsi 

n^ m \ m m> J 



m* 



m — I ^ 



(5o) { ^is^s— ^8*+i 

2 



2 \ 3 a 2 / 






Nous avons, d'autre part, 

(5i) i-h8=(i + Ay 

d'où 

o=(i-hAr-i= - A-+-f-^ AÎ-+-... 



et aussi 



par conséquent, en remplaçant m par — dans (5o), 

(52) m*8^ = A^-h^m/^ — j A^+i-+- T 5m( — j -i-5,m«/^— V 1 A*-^*-4-. . . 

ou bien 

^ A* m — I A*-»-* m — \ (im — i s — i m — 1\ A*+* 

8* = s 7-7 4- s 



\ (im — I s — I m — 1\ 

\ 3 2 2/ 



m* 2 m*+* 2 \ 3 2 2 / m*-^-* 

Si A'' est la dernière différence sensible (A'^*==o), on aura 
évidemment 



o'' = — A'', 
m*' 

8r-l = \^r-\ _ i^ L 1 Ar 

m'-» 2 m'' 



\) l'i.m — I r — 3 m — 1\. 



(r — 2)(/?? — \) l7.m — I r — 3 m — ^ \ * ;. 
2 m'' 
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En faisant s=: i, 2, 3, ... dans (52), on trouve 



= — A A* 

m 2 m* 

(m — \)(Sm — 4) ^3 ^ (m — i)(3m — i)(4m--2) ^^ 

m* 2 /n* 

2(m — ï)(ii/n — 7) 2(m — 0(4'^ — 2)(5m — 3) 

>, 1 ., 3(771 — 1) . • 

o3 = — A* ^ A* 

3(m — i)(i4m — 10) ,, 3(m — i)(5m — 3)(6m — 4) ,. 
24 m5 48 me 



Le cinquième terme de S est + — — — -^^— ^ — LM^! — 12 A^. 
^ 120 m* 

L'expression des différences 8 que C.-F. Degen a donnée sans 
démonstration dans une courte Noie (Astr. Naeh., n° 15, 1822), 
cesse d'être exacte pour 5 > 3. Le Verrier a traité le mêm^ sujet 
(Ann. de VObs,, t. I, p. 126 et i5i). 

H. Les différences S sont les différences successives à partir 
de fo(8^oj S^fjj, . . .), comme les A sont les différences successives 
à partir de To, et elles servent à former, par additions répétées, 
les termes cherchés t^y t^^ . . •• Dans le Tableau, ces différences 
descendent en diagonale. Lorsqu'on veut faire usage des différences 
médianes D, il faut construire un autre algorithme. 

Remarquons d'abord qu'en appliquant aux ? les formules (4^^) 
il vient 

i 2é; = i + 5-(i + S)-i, 
^ I c?2 = i-+-8 + (i-4-8)-i — 2, 

où rf, rf^, ... sont les différences médianes pour les intervalles 

réduits. Or, en vertu de (5i)^ nous avons 1 4- S = (i + A)'"; et, 
en faisant [jl === — > (48) nous donne 



m 
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En substituant dans (54)) on trouve 






î7l-»-l 

(56) 






V 

OU bien 



(57) 



'^ V/^ /3 V^^ /S 

= -L D - ^'""' D3 + (^'-0(4^^-1) Ds 
m 6m3 i2o/w5 



et 



(58) 






De ces expressions fondamentales, les autres d se déduisent par 
les relations 

(69) < 

d.»-. = rf.rf.» = —1- D««+« - ("■+^K^' — ) D,„^, + . . . . 

On aurait, par exen^ple, 

7»^ 4 m* 120m' 

7?i» 6 m^ 240 /n* 

Il ne faut pas oublier qu'on n'a pas d.d^=d^f mais qu'on a. 
d'Kd^ = d'^-^K 

Ces formules sont d'accord avec celles que M. J.rJ. Astrancf'a 
publiées en i858(yi. N., 1125) et en 1875 {Vierieljahrschrift, 
t. X, p. a^g), m.ais sans en indiquer complètement la loi. Encke 
avait traité un problème analogue dans le Jahrbuch pour i852. 

En suivant une marche inverse> on pourrait exprimer les D par 
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les rf, et Ton trouverait 

D = mû?-h(/îH-i)8rf'-*-(/w-+-a)5û?«-H..., 
Pour /n = a, on aurait donc 

Ces relations peuvent faciliter le calcul des difTérences d. Ayant 
calculé ces dernières pour l'argument o, on trouve les différences 8 
par les formules 

i.,=d-\d>, !rïi;~^'" 8î,=d»-ws 

«.=''-iA 8r=>v8i. '-^-^'^'^^'- 

J'ai dû me borner à démontrer ces diverses formules en renon- 
çant à donner des exemples de calcul, qui auraient pris trop de 

place. 

Quand nous ferons d'autres applications des différences mé- 
dianes, nous désigfïérons par la lettre M la moyenne ^ ( A H j , 

représentée plus haut par D, et D représentera l'opération fic- 
tive ' On aura M = D y/i-f-ÏD^, et, pour unifier la notation 

des différences successives, il suffira de diviser par y/i -f- ^-D^ les 
développements où figurent D2«+*(o) ou D2'»(i). 

12. Avant d'aller plus loin, il convient de faire une remarque 
qui peut paraître banale, mais qui n'est pas superflue. C'est qu'au 
lieu de la quantité y dont on connaît une série de valeurs (soit 
isolées, soit relevées sur un tracé graphique), on peut interpoler 
telle fonction de j^ qui s'y prête avec plus de facilité. •On peut 
essayer j^*, ou ^y^ ou logy. Dans beaucoup de cas, surtout dans 
les problèmes de Physique ou de Statistique, c'est le logarithme 
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t. 84, 1877), M. Hermite indique une autre généralisation de la 
formule de Lagrange, qu'il obtient en cherchant à représenter par 
un polynôme F du degré n — i, non seulement la fonction 
donnée f{x)^ mais encore ses premières dérivées, de sorte que, 
par exemple, /(a), /'(«)? •••, /*"'(<^) coïncident avec F(a), 
F'(a), . . ., F«~*(a). On peut consulter, à ce sujet, une Note de 
M. V. Williot {BulL des Sciences math., octobre 1890). 

Pour étendre la formule de Lagrange au cas de deux variables ;r, 
y^ on sait qu'il suffit d'assimiler y à x^ en posant 

où les facteurs X/ et Y/ sont de forme semblable. 

Il y a quelques années, à propos d'un travail de M. Weyer 
(Astr. Nackr., 2804), j'ai montré comment la formule de La- 
grange se transforme quand le nombre des ordonnées yn devient 
infini, et en particulier quand ces ordonnées forment une suite 
périodique où un groupe de m ordonnées se reproduit indéfini- 
ment; les coefficients se changent alors en fonctions trigonomé- 
triques (voir BulL astr,, t. IV, p. 5i5). 

Mais je laisserai de côté l'interpolation par les séries trigono- 
métriques, dont Lagrange, on le sait, s'est occupé à plusieurs 
reprises (Recherches sur la manière de former des Tables des 
planètes diaprés les seules observations, 1772; Œuvres, t. VI, 
p. 507; Sur les interpolations, 1780; Œuvres, t. VII, p. 54 1). 
11 y aurait à citer la formule de Walbeck (Corresp. astr., 1820), 
celle de Bessel (Astr, Nac h r,, n"^ 136; 1828); puis la méthode que 
Le Verrier a développée dans le tome I des annales de V Obser- 
vatoire (1855) et dont, peu de temps après, Encke a donné un 
exposé dans le Jahrbuch ( 1860). M. Hoiiel s'en est occupé à son 
\.owv (Comptes rendus, 1861; Ann. de VObs., t. VIII, i865). 

Le problème a été résolu par Gauss, d'une manière beaucoup 
plus simple, dans le Mémoire posthume intitulé : Theoria inter- 
polationis methodo nova tractata (Œuvres, t. III, p. 265), où 
l'on trouve la formule 

. X — b . X — c . X — l 
sin sin • • -sin 

'X 1 '1 

^ . a — o . a—c . a — l 
sin sin . • -sin 

'2 X '}. 
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qui peut tenir lieu de la série 

y = oi-\- ol'cosx -h oi^cos^x -4- 
P'sina? ■+- ^"81020? 



II. 



13. Dans beaucoup de cas, la recherche des coefficients d'une 
formule d'interpolation se ramène à la résolution d'un système 
d'équations linéaires. 

L'introduction des déterminants est un excellent moyen d'élu- 
cider la théorie de l'élimination; mais, lorsqu'il s'agit de calculs 
numériques, on préférera toujours le procédé ordinaire de l'élimi- 
nation successive des inconnues. Ce procédé peut d'ailleurs être 
présenté sous une forme élégante, en suivant les prescriptions de 
Gauss ou celles de Cauchy. 

Supposons d'abord qu'on ait à déterminer m inconnues or, y, 
Zj . . • , ^ par m équations linéaires 

Aia7H-Bi^-l- Ci-5 +.. .-4- H|«= Fi, 



dont nous représenterons l'une quelconque par l'équation typique 

(I) Xx-hBy-\-Cz-i-,.,-\-Ut = F. 

On commence par éliminer x, en retranchant de chaque expres- 
sion F la première F| , multipliée par le rapport ^ des coefficients 
respectifs de x, et l'on obtient m — i équations de la forme 

(II) j^AB-h^AG4-...-h«AH== AF, 

les diiTérences premières A étant définies par les relations 

AA=o, AB = B-Abi, aG = G-Ac„ ..., aF = F-Af„ 

Al Al Al 

de sorte que ABi, ACi, . . . , AFj s'annulent évidemment. L'élimi- 
nation dey entre les expressions AF2; AF3, . . • , AFm donne m — 2 
équations nouvelles du type 



— 26 — 
en faisant 

A«B=o, A2G = AG-^AGî, ..., AîF = AF-^AFj, 

Ai32 Adj 

et ainsi de suite. Il est visible que AF| ; A^Fi, A2F2; .. , 
A^'Fi, . . ., A^'Fto sont nulles par définition. 
On arrive finalement au système 

S' xA.-hyB-[-zC-^,..-^ tll =F, 

7 AB -j- ^ AG H- . . . H- ^ AH = AF, 



t A«»-i H = A'«-» F, 

dont la première ligne représente m équations, la seconde m — i, 
et la dernière une seule, comme le montre le Tableau suivant des 
seconds membres 

/ Fi, Fj, F3, . . ., F„i, 
AFj, AFj, . . . , AF,;i, 

i'^) { A2F3, ..., A* F;;,, 



. • • • . 



A'«-i F 



/;i- 



Pour le calcul des inconnues x^y, 5, il suffit de prendre la pre- 
mière équation de chaque type, à savoir celles dont les seconds 
membres soilt 

(3) Fi, AF2, A2F3, ..., A'«-iF,„, 

en remontant de là dernière à celles qui la précèdent; elles four- 
nissent ainsi, l'une après l'autre, les valeurs de ^, . . . , >s, y, x. 

Si les inconnues sont liées entre elles par n équations de la 
forme 

( «1 a? -h 61 j^ -+- Cl >5 -h . . . -h Al / = /i , 

(4) < aiX-h bijr-h CiZ -h,..-^ hit =zfi, 

f , 

dont le nombre surpasse celui des inconnues (n > m), ces équa- 
tions sont, en général, simplement approchées, et ne deviendraient 
exactes qu'en remplaçant f par f — e, où s représente l'erreur 
dey*. Il faut alors commencer par en déduire un système d'équa- 
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lions finales en nombre égal à celui des inconnues; ce qui peut se 
faire de deux manières. 

On peut, d'une part, en ajoutant les équations (4) après les 
avoir multipliées par des facteurs )vi, ..., \nt ^u •••> J^/ii •••? 
former un système de m équations linéaires tel que le système (i), 
et en déduire les équations finales (3). En partant de Téquation 

type 

ax -^ h y -\- c z -h . . .-h ht = f, 

et nous servant de la notation de Gauss 

(ap ) = ai p, -f- aj P, -^ . . . 4- «;, P„, 
les sommes pondérées des équations primitives pourront s'écrire 

('ka)x -^('kb)y -\-{'kc)z -|-... = (X/) = F,, 

(vrt)a7 -^-(yb)y -^(vc)^ -i-... = (v/) = F3, 



L'élimination successive de ar, de y, de z, ... donne des séries 
d'équations telles que les suivantes 

( |i6. 1)7 4- f [xc. i)>s -f-. . . = (|x/. 1), 

. ; (vÔ.l)j^ -h(vC.l)3 4-...= (v/.l), 

(VC.2)3 -+-. . . = (v/.2), 



• ■ • • 



où, d'après la notation de Gauss, 

(>/.,) = (Hi/)-i|Jj(X/), ... 

Avec la notation de Cauchy, qui représente par A, A*, . . . des 
« espèces de différences Gnies de divers ordres », ces expressions 
deviennent 

AF, = A((x/)=(,V)-([^a,)j4^' 
A» F, = y-(',/) = A(v/) - A(v6) ^^J, 
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et il est visible qu^ellcs peuvent être mises sous la forme 

AF, = (fiA/), A«F8 = (vA«/), 



• • • > 



en posant 



« 
ce qui, à la vérité, modifie légèrement le sens du signe d'opéra- 
tion A. 

Il faut se garder de confondre le symbole A, employé ici, avec 
celui des différences finies ordinaires, où A. A = A^. Les A sont ici 
des coefficients d'élimination. Ainsi A^ naît de l'élimination de x^ 
y entre les équations 

F = Aa? -f-Bj -h..., 

Fi = Aiar-i-Bij'-h..., 

Fj = ktx 4- Bj^ H-. . . , 

et l'opération se fait comme si les autres inconnues n'existaient 

pas. En désignant par [a6|] le déterminant formé avec les élé- 

a b 
ments , on trouve 

AF =F-^F.= [^, 
Al Al 

AB, '^^ [AiB,]' 

A3F = A.F-^î£A«F3=[i:ALÇl^, 

A«G3 ' [AiBjGa]' 

et ainsi de suite. On voit que les A sont des fonctions symétriques 
des coefficients qui y sont contenus. 

Les quantités ((x Ay*), (v A^/), . . . sont composées d*une manière 
analogue, les éléments des déterminants étant les coeffîcients()^y*), 

Il nous semble que l'emploi du signe A aide à régler la marche 
des opérations. 

On se trouve ainsi conduit (comme l'a remarqué Bienaymé) au 
mode d'élimination fondé sur l'algorithme de Cauchy. 

14. On peut, en effet, avec Cauchy, établir directement les 
équations finales sans passer par toute la série des. expressions (5). 
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Pour cela, remarquons d'abord qu'en éliminant successivement x^ 
y, Zy • . • des équations (4), on peut écrire à la fois le système des 
relations approchées (où chaque ligne représente n équations) 

xa-^-yb -^ zc -4-...-4-f/i =/, 

yi^b-^-zS^c 4-...-hfAA = A/, 
(7) { ^ A2c -+-... -4-«A« A =A2/, 



t A"»-» h = A'»-»/, 



et le système d'équations définitives 



xCka)-hy{lb) -+-^(Xc) -h.,,-ht(kh) =(X/), 
, . 7(|iA6)-4-^({i.Ac) -4-...-hf({iAA) =(fiA/), 

^(v A2c)H-. . .-f- *(v A* A) = (v A»/), 



qui se déduisent l'un de l'autre, ligne après ligne, les différences A 
étant formées comme il suit : 

v=/-[Mi-. ^-/-v-JJIg" 

La première des équations (8) est la somme pondérée des n 
équations contenues dans la première ligne de (7); combinée avec 
celles-ci, elle sert à former les n équations que représente la 
seconde ligne de (7), et dont la somme pondérée fournit la 
seconde équation du système (8). Celle-ci, à son tour, sert à 
déduire la troisième ligne de (7) de la seconde 5 cette troisième 
ligne fournit la troisième équation (8), et ainsi de suite. 

Pour abréger, nous introduirons les notations nouvelles ( * ) 

„_(X/) >_ (f^A/) ,^ (vAV) 

(9) \ ^"(^«)' ^ (f^A*)' ^ (^A*0' 

^f = f--^a, A«/=A/-cp'A^ A3/=A«/-?''A2c, 



(*) Il faut dire que la notation que nous avons adoptée n'est pas tout à fait 
celle de Cauchy, qui écrit, par exemple, 

A/ = /-a(V), a= "" 



mais les formules gagnent en symétrie par l'introduction des rapports 9. 
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qui entraînent ces formules à* opérations 

(XA) =0, (XA«) = o, ([iA2) = o, 
(lo) {(XA8) = o, (fjiA8) = o, (vA3)=o, 



On a, en effet, par définition, 

(XA/) = (X/)-cp(Xa) = o, .... 

En désignant toujours par a, p, y, . . • les rapports o pour 
f=. a^ bj Cj . . . , on trouve 



a = 1, 



(II) 



__ (Xc) ,_ ([xAc) ,__ 

^ "(Xa)' ^ ""([xA^)' ^ -' 



(12) 



Aa =0, 

A6=6 — Pa, A2^»=o, 

Ac = c — Y^> A^c = Ac — Y^^i A3c = o, 



le système (8) peut s'écrire 

z -I-.. .-+- <t''= cp", 



et les formules (9) donnent immédiatement 

1/ = cpa -H cp' A6 H- cp" A^c -I- . . . , 
A/= cp'Aô -H cp''A2c -+-..., 

A2/= cp"A2c-h..., 

15. Le procédé d'élimination qui vient d'être exposé donne les 
équations finales pour la détermination des inconnues sous une 
forme très simple. Ainsi que l'a fait déjà remarquer M. Bienaymé, 
il convient également à la méthode des moindres carrés et à la 
méthode de Cauchy : la différence est dans le choix des facteurs 
X, [JL, V, ... qui servent à former les sommes pondérées des équa- 
tions approximatives (7). 
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On commencera par supposer 

.en faisant^' = o, x; = o, . • . , et le système (i3) donnera 

m ' 

en première approximation. Dans la méthode de Cauchy, le rap- 

port <p = W^ devient 

S/ 

9= c^' 
^ Sa 

où Sa=.^±a est la somme des valeurs numériques de a, et 
S/= S ±/, les changements de signe du numérateur étant réglés 
sur ceux du dénominateur, dont tous les termes deviennent posi- 
tifs. Ensuite on formera les résidus 

dont la combinaison (XA/) ou SA/ doit s'annuler, ce qui fournit 
une vérification du calcul. S'il y a lieu, on procède à la seconde 
approximation, en posant 

et le système (i3) donnera 

r = ?'. 

OÙ 

Les relations (i4) donnent d'ailleurs 

On peut continuer de cette manière, jusqu'à ce que les résidus A'"/ 
paraissent suffisamment petits. 

17. Présentée de cette manière, la méthode de Cauchy se rattache 
très naturellement aux procédés d'élimination plus anciens, dont 
elle constitue une ingénieuse modification. Je n'ai pas voulu en 
interrompre l'exposé par des remarques relatives à l'historique de 
la question, qu'il y a lieu de placer ici. 

R. 3 
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tuées, et de remplacer «, i', (v, . . . par 

/yO >/.' ^" r« 

Dans le cas le plus simple, celui d'ua développement suivant les 
puissances ascendantes de x, on aura donc 

/y»0 — 1 'T»' — />• /y»" — 'y*? 

et :r'* représentera, à l'ordinaire, la puissance n^"^^^ de ^; mais, en 
général, nous désignerons par a:" la fonction de x dont n est le 
numéro d'ordre ou le rang dans la série. On verra qu'il n'en peut 
résulter aucune confusion et qu'on peut se dispenser d'écrire x^"'^ 
au lieu de ^", ce qui serait peut-être plus clair, mais très incom- 
mode pour la transcription des formules. Nous conserverons aussi 
la notation de Gauss 

pour représenter les sommes de produits ou sommes pondérées. 

19. Ceci posé, nous allons considérer le cas spécial où le pro- 
cédé de Cauchj est appliqué à la construction d'une formule d'in- 
terpolation. 

En admettant que la fonction de x accessible à l'observation (*) 
soit développable en série convergente suivant les puissances de x 
ou, plus généralement, suivant certaines fonctions x^^ x'^ x"j . . . , 
de sorte qu'on puisse poser 

(i5) /= Aa70-^B:r'-4-Ca7"-h..., 

il s'agit de déterminer les coefficients A, B, C, . . . à l'aide d'un 
certain nombre de valeurs données /i, f'^^fzy . . ., qui correspon- 
dent aux valeurs x^^ x^, X2, ... de l'argument x. On aura donc 
le système d'équations linéaires 

fi = Xx\ -\-Bx\-\-Cx'[-h,.., 
/2 = Xxl -h BiPj ■+- Gxl-\-.. ., 



(') Dans beaucoup de cas, il sera plus commode de développer directement 
quelque fonction de /, par exemple son logarithme, comme nous l'avons déjà 
dit (n» 12). 
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oii les coefficients donnés xf ^ x'^ ^ x'I j . . . sont ceux que nous 
avons désignés plus haut par a/, 6/, c/, . . . , tandis que A, B, C, ... 
remplacent les inconnues x^ y, z, .... 

Gomme le nombre des équations est, le plus souvent, supérieur 
à celui des inconnues, on déduit du système donné un nombre 
d'équations strictement suffisant en formant les sommes pon- 
dérées 

(i6) * (>/)=(X^o)A-+-(Xa7')B-t-(X^^)CH-... 

analogues à celles qui composent le système (5). 

Mais le nombre de ces équations dépend évidemment du nombre 
de termes qu'on veut conserver dans la formule d'interpolation (i 5). 
Or, on sait rarement d'avance combien de termes il faudra intro- 
duire pour obtenir une représentation satisfaisante, et pour chaque 
terme nouveau on est obligé de répéter tous les calculs^ C'est 
pour parer à cet inconvénient que Cauchj a imaginé sa méthode 
d'approximation successive, où lés termes nouveaux se déter- 
minent en gardant toujours les anciens. Pour y accommoder nos 
formules, nous n'avons qu'à mettre x^^ x' ^ x" ^ ... à la place de a, 
^, c, . . . , et A, B, C, ... au lieu àe x^ y^ Zy . . . , Les relations ( 1 4) 
montrent que/* s'obtient sous la forme 

(17) /= cpa70-i-cp'Aa7'-f-cp''A2a7''-+-. . ., 

les coefficients co et les A étant définis par les relations suivantes 

(18) A/ = /-cp^o/ AV = A/-o'l^', A3/=:AV-?''A*^", ..., 

^'^^ ^-(I^)' ^-ÔIÂ^' ^-(TÂi^F)' •••' 
qui entraînent les formules d'opération générales 

( (XA) = o, (XA2) = o, (hlA2) = o; 

(10) 

f(XA3) = o, ({jiA3) = o, (vA3) = o; 

déjà indiquées plus haut. Ajoutons que, n étant le nombre des 
valeurs données de/, on aura A"/=o. 

Les rapports cp, qui étaient dénotés plus haut par a, P, y, ... 
pour /= 67, ^, c, ..., seront maintenant désignés par (po, 'fi, 
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Sj, ... pour/^ x", x', x" , .... On aura donc l'algorithme 



Aa:<> = o, 

Aa?' =3?' — cpi^P^, A* 47' = o, 

(21) ^' Aar' = ^' — cpjaro, A'^c" = Aa?" — îp^ Aar', A^a:" = o, 

Aa7"'= ar'"— cpaaro, A«;r"'= Aa?'"— cp'3 Aa?', A»ar'"= A«ar"'— cpj A2ar% 



a?o= a?», 

, . F a;' = ©laro-h Aar', 

(22) { ^ 

x" = cp,aro -h cp'g Aar'-f- A^a?', 



On commence par poser y*=A:r® en première approximation, 

et Ton trouve 

A = cp, /=çpa?o. 

On calcule alors les résidus A/*, à savoir 

et s'il y a lieu, on procède à la seconde approximation en posant 

/= Aa-o+Ba?'. 
La formule (17) donne alors 

/= ça^o-f- cp'Aa-', A/= cp'Aa:', 

par conséquent 

B = cp', A = cp — cp,cp', 

OÙ les coefficients cp sont' fournis par les formules (19) et (20), que 
nous pourrons écrire, avec Gauchy, 

_ Sa-^ S/ ,_ S^A/ 

'^*~ Sa^o' "?" Sa^' "^ " S'Aa-'' 

en réduisant les facteurs )v, [jl à =h i , les signes étant pris de ma- 
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mère que tous les termes de S:r*^ et de S'A^' deviennent positifs. 
On dira alors, avec Villarceau, que x^ est la dominante des 
sommes S, et Aj;' celle des sommes S', ces deux fonctions détermi- 
nant, l'une le facteur a dans les sommes S, l'autre le facteur [jl 
dans les sommes S'. 

On cherche ensuite les résidus 

A«/= A/— ç'Aa:' 

et, s'ils sont encore sensibles, on essaye une troisième approxima- 
tion 

dans laquelle on a 
et, par suite, 

les coefficients étant toujours donnés par les formules (19) et (20). 
On peut continuer ainsi, jusqu'à ce qu'on arrive à des résidus A'"/* 
qui soient négligeables. Mais il est visible que le choix des fac- 
teurs X, [JL, V, ... reste arbitraire. 

Au lieu de les réduire, avec Cauchj, à dz i, nous pouvons les 
égaler aux coefficients des inconnues A, B, G, ... et prendre 
\ = x^j [JL = j?', V = x", . . . , pour nous conformer à la méthode 
des moindres carrés et obtenir le résultat le plus avantageux. 

Il est toutefois permis de penser que l'emploi rigoureux de la 
méthode des moindres carrés est bien rarement nécessaire, et 
qu'il suffit, le plus souvent, d'en suivre le principe en arrondissant 
les facteurs, c'est-à-dire en remplaçant des groupes de facteurs 
par une valeur moyenne en nombre rond. Cauchy, en les rédui- 
sant tous à ±: I, va probablement trop loin; on sera plus près de 
la règle, ainsi que nous l'avons déjà dit, en égalant les facteurs 
ào, dz I, ±: 2, ... selon l'importance du coefficient de l'inconnue 
qui joue le rôle de dominante, 

î20. Dans le cas le plus simple, celui d'un développement sui- 
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^ant les puissances de Xy nous avons 



X ^— I f X — 37. X — ^ X f • • • ^ 



et (17) devient 



/= A-i-Ba?-4-Car«-t-... 



<23) 



y= <p -+-çp'AiC-+- «p'A'aj^-f- 



Les facteurs X se réduisent à +1 (dans la méthode de Gauchy, 
comme dans la méthode des moindres carrés); la somme S est la 
somme ordinaire S, et (X:c®) = n représente le nombre des obser- 
vations; les quantités cp, cpi, ©2, • • • sont des moyennes arithmé- 
tique^. On a 



(24) 



?3 



= ^V, 



n 



1 



n 
n 



- 2ar3, 
n 






?l = '» 



?1 = 



<?»= 



( (JL A^* ) 
([JlA3?3) 

(fxAa:) ' 



(vAV) 






T3 



(V A23?3) 

(vA«:i-«)' 



• » 



„/// , 



Ensui te 



>5> 



Ao? 

Aa:* 

A37' 



= r — co, 



/ 

a?2 
X^ 



?1» 
?3, 



A* a? 

A2^2 

A2a73 



= A/ - cp'Aa:, A3/ = A*/ - cp'A'o;», 

= 0, 

= Aa?^ — (p2 Aa7, A3 3?* = o, 

= Aa:3 — cp3 Aa?, A3 373 = A' 073 — ç^'^ A*^', 



• • » 



A*ar3 = 






enfin 



(9.6) 



I =cpo, 




37 — ©1 -i- A37, 




372 _ O2 -h O2 A37 -H A2372, 




a73-_cp3-+-'^'3A37-f-cp'iA2 372H- 


A3r^ 
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La forme explicite de la série devient 

/= '^ -f- çp'(a: — çp, ) -f- cp"[ a:2 _ cpî — cp; (a; — ©1 )] + ... . 

Quant aux facteurs |jl, v, . . . , on prendrait, avec Cauchy, 

(v A2a72) = S*A«a?2 = 2 ifc A«iF*, 



en choisissant les signes de manière à n'avoir ici que des termes 
positifs. C'est sous cette forme que Bartlett a étudié le procédé 
d'interpolation de Cauchy, en traitant quelques exemples numé- 
riques (^m. J. of Se, 1862). 

Si nous voulions nous conformer au principe des moindres 
carrés, il faudrait prendre [jl = ^, v = x^^ .... 

21. De même, dans le cas général, en faisant usage de la for- 
mule 

(17) /= (pa:0-t- cp'A^'+ fl^x"-^. . . , 

on se conformerait au principe des moindres carrés en prenant 
X = :r®, [ji = j:', v = x"^ . . . , et l'algorithme aurait la forme sui- 
vante : 

^f = f-^x\ A2/=A/-cp'Aa7', 

AV=A2/-o''A2:p", ..., AV=o; 

^ (x^x^y ^ {x'^x'y ^ (x"^^x"y 

On aurait, pour/= x^^ 
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par conséquent 



(aro^'^) , (x'ix") 



La formule 



donne d'ailleurs 






en déHnissant cp(*) par la relation (a;* A*+*y) = o, ou (x^~* A*/) = o, 
c[ui entraîne (:r^"* A^"*'*y) = 0, . . ., et enfin, d'une manière géné- 
xale, 

<28) (x'^.\^) = o pour r < 5, 

comme le demandent les formules (10), qui s'écrivent maintenant 

(10) (J70A) = 0, (a70A2)=o, (a7'A2) = o, 

On a même, plus généralement, 

(29) (A"*47'*.A*) = o pour r < 5. 
Cela résulte des égalités 

(30) (F.A*/) = (AF.AV) = (A«F.AV) = ... = (A^F.AV), 

qui se déduisent successivement des formules (27) et (28); elles 
montrent que, dans les formules (28), on peut remplacer x^ par 
^x''j A^x'*, . . ., à^'x'' et en général par Ùl^x^^ puisque cette dilTé- 
rence s'annule pour m >> /*. On aura donc aussi 

(3i) (\''x'',l^xs) = o pour r^s. 

Grâce à ces relations, les sommes qui figurent dans les expressions 
des coefficients cp peuvent s'écrire de diverses manières, par 



exemple 



(32) 
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{x'\f) =(flx') =(A:p'.A/), 

(x'^x') ={^x'.^x'), 

(x" \^x") = i^xW^x") = (l^x'.ùi^x'' ), 



et l'on voit que les dénominateurs des cp deviennent des sommes 
de carrés. M. Carvallo a fait une remarque analogue. 

' 22. Considérons en particulier le cas d'un développement sui- 
vant les puissances de x. Nous avons alors x^ = i, x'= x^ ... et 
(x^x^)= n. L'algorithme, pour l'application de la méthode des 
moindres carrés, devient 



/ / =0 + 0' 



A^c + cp'A'iT*-!-.. ., 



(33) IV = o'Aar-hcp'A^arïH-..., 

( Aî/= fl^x^-h.,., 



pUl! 



(34) 



A/ =/ -cp, Aî/ =A/ -^'\x, A3/ =A2/-cp''A«a:î, 

\x =x — cpi, A2.r =0, 

Aa72 = a7* — cpî, A2a7' = Aa?* — ^t^^? A3a:*=o. 



et 



n •^' * (:fA^) • (a72Aîa72) 



T» = T7-^» cp; = i, 



/i 



(35) { 1 , , , (xiix^) 

1 , (x^x^) „ (37»A»ar3) 

<P3=-1^3, Ç3=__, ,p3=___, cp3 = I, 



et aussi 



a: = çpi -h Aar, 

2g ) 372 = ©2 -f- ç>'^ >^a7 -h A2a:-2, 

i .r3 = Cp3 -h Cp3 Aa7 + Cp'i A2j-2 4. A3;r'\ 
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11 est à remarquer qu'en vertu des relations (32) nous avons 

{x\f) =(/A:F) = (Aa7.A/), 
(a?A2?) =(Aa7.Aa7), 
^, j (rpAo?») =(x^,^x) =(Aa7.Aa72), 

' (ar2A2/) =(/.A2a:2) = (A2a?2.Aî/), 
(a:2A2a72)=:(A2ip2.A2a^2), 



ce qui permet de transformer les coefficients cp^-^^ de bien des 

iTfianières. 

Si nous voulions les exprimer par les sommes de puissances 
^i7 'fa? ?3î • • •? il suffirait d'appliquer aux relations (36) les for- 
mules (28), qui peuvent s'écrire 

(38) 1:A=o, 2A» = o, 2a7A2 = o, ..., SiF'-A''-^'' = o. 

On trouverait, par exemple, 

/i(cp2— cpi9i) = (a?Aa7), 7i(çp3 — cp20i) = (a72Aa7), 

7i((p3 — Oicp2) = (a:Aa7)cp;, Il{t!^,,— o^^^) = {x^ ùiX)^'^-\- {x^ ^'ix^), 

n(cp4 — o,cp3) = (a7Aar)cp;,, /i(cp5— cp2?3) = (^* Ajc)o'3 -{-(a:» A2a72)çpJ, 

et ainsi de suite. 

23. On pourrait aussi exprimer les coefficients delà formule (33) 
par les puissances de A^. Posons 

Aa? = >s, a7 = 5-l-cpi, 

nous aurons, en ve/tu de (28), 

(89) 1:^=0, (-3*-'A0 = o, (5*A^)=(ir^AM, 

et ces formules permettent de déterminer les coefficients a,, p,, 
y^, ... de l'expression 

en fonctions des sommes S:;"', par une série de relations de la 
forme 

(iO j 0= Ss^+i + a^Si-'^ -h p^S^-'f-ï-h Y-^^-^"' 



r.Ç-î 
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Elles donnent encore 

et 

Or la différence z = x Hx ne change pas, si nous rempla- 
çons X par :c + «, et A/* ne change pas non plus, si l'on mety-f-è 
à la place de /. Il s'ensuit [comme le montrent les relations (34) 
et (3^)] qu'il est permis, dans la formule (33), d'ajouter une con- 
stante quelconque à x aussi bien qu'à/; on peut aussi remarquer 
qu'elle est homogène (de dimension zéro) par rapport à x, de 
sorte que Tunité de x reste arbitraire. 

Quand les valeurs données fi y /2» ••• correspondent à des 
intervalles égaux de l'argument x ou seulement à des valeurs de x 
distribuées s^^métriquement autour de la valeur moyenne, on aura 
Zi = — Zfi^ et en général zi= — Zn^i_i', par suite 

(42) 23=0, 2^3 = 0, 2^5 = 0, 

et les équations de la forme 

0= a^2i2'w_|_ Y5 252/«-2_|_£^2^2m-4_j_._ 

qui résultent de (4i) donnent alors 

«5 = 0» "(s = 0, £5=0, . . . , 
de sorte qu'on trouve 

(43) ^'x^ = z-'-h ^s^'-^ -H- ô^^^-* + . . . . 

Gomme l'unité de x reste arbitraire et que les formules ne ren- 
ferment que des différences, nous sommes libres de supposer que 
Xi^ X2, . . . représentent les nombres i, 2, 3, . . . , /i, de sorte que 

,,,^ n-4-i n-4-i . n-\-i 

(44) ^1 = > z = x 9 Zi=i 

* a 9. 2 

et 

^^ n — I n — 3 

Zi = iiz > IlZ > .... 

2 2 

On a aussi, d'une part, 

;r* = (©!-+- s)^^ = o'j ■+-. . .H- 5cpi5^-iH- 5-% 
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et de l'autre, par les formules (36) 

et, en tenant compte de (43), la comparaison donne 

1 ?^ = 1.2 ??-P^> 



24. La formule (33), à laquelle conduit Talgorithme de Cauchy, 
accommodé au principe des moindres carrés, ne diffère que par les 
notations de celle que M. Tchébychef a donnée en i854 et dont 
on trouve la démonstration dans son Mémoire de i855 sur les 
jf raclions continues^ traduit en français par M. J. Bienajmé, 
en i858 {Journal de Liouville, 2® série, t. III). Cette formule 
peut s'écrire 

f=^o(^)-— H^i(a^)— 1-..., 

où les fonctions ^s sont celles que nous avons désignées par A^^^, 
tandis que 9^ représente le poids des observations, que nous avons 
pris égal à l'unité. 

M. Eug. Rouché a étudié les séries de ce genre dans un Mémoire 
sur le développement des fonctions (Journ. de V Ecole Polytech- 
nique, i858). M. Hcrmite, dans une Note très intéressante sur 
V Interpolation {Comptes rendus, 10 janvier 1859), rattache les 
résultats de M. Tchébychef d'une manière directe à la formule 
de Lagrange, généralisée par l'introduction des polynômes 8(^). 
M. Borchardt s'en est occupé, à son tour, en 1860. M. Tchébychef 
lui-même y est revenu à plusieurs reprises : d'abord dans le Mé- 
moire sur V Interpolation par la méthode des moindres carrés 
{Mém. de VAc, de Saint-Pétersbourg , 1869), ensuite dans deux 
Mémoires sur l^ Interpolation, publiés en 1864 et en 18^75 (en 
russe, dans les Bulletins de l'Académie), modifiant chaque fois les 
notations et le mode d'exposition. Il a aussi traité le cas où la 
fonction y* est en partie connue, de sorte qu'il s'agit de la repré- 
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î25. La formule de M. Tchébychef, destinée à représenter des 
observations de même poids, lorsque /est développable en série 
suivant les puissances de a:, peut s'écrire comme il suit : 

(46) f=u^)—- — -^i(^)== — -.... 

Les fonctions ^s{^) sont les dénominateurs des réduites qui se 
rencontrent dans le développement en fraction continue de la 
somme 



(4/) 



ju ** 1 ^ «^2 "^ ■ '*^ l\ 



On a d'ailleurs ^j;o= i, et l'on s'assure facilement que la fonction 
^s{^x) de M. Tchébychef n'est autre chose que la différence A* a? 
de Cauchy. En vertu de (3i), on a donc 

(48) 24;;..^;,==o. 

La formule (46) revient à celle-ci 

qui coïncide avec notre série (33). Mais M. Tchébychef y arrive 
par une autre voie, et fait usage d'algorithmes ingénieux pour le 
calcul des coefficients. Dans le Mémoire de 1869, leg sommes 
(y*. A*^*) sont ramenées à la forme équivalente (^*.A*/) et les 
calculs reposent sur des séries récurrentes; on les applique à un 
exemple numérique. 

26. Un cas particulier très intéressant est celui où les valeurs 
données yi correspondent à des intervalles égaux de l'argument x. 
Comme nous l'avons déjà dit, il est alors permis d'idenlifier Xi, 
^2> ^:m • • m «^m sivec les nombres i, 2, 3, . . ., n, de sorte que la 
série (47) devient 

(49) 



X — I X — 2 X — Il 

K. ', 
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En désignant par D* les différences finies ordinaires, M. Tché- 
byclief montre qu'on aura 

(5o) ^s(^) = \y'^sy 

où 

(5i) 4>5 = (a7 — i)(a? — i),. .(x — s)(x — n — i).,.(x—n — s), 

à un facteur constant près, qui reste évidemment arbitraire. 

La fonction ^s, ainsi définie, joue dans le calcul inverse des dif- 
férences un rôle analogue à celui des polynômes de Legendre dans 
le Calcul intégral. 

La formule (46) se transforme comme il suit : 

■ 

(5.){ • «;. . 

5 , , .'^ i(i-h\)(n — iXn, — i — i)^^. 

1 

Dans le Mémoire de 1809, M. Tchébychef remplace œ par z = A.r; 

on a 

n H-i 

X = z-\ , 

2 

par conséquent 






1S 

'À 



ou 



bien 



'^^-»\/^.__3.-"'-^^ 



4 



/i2 — I 



Comme on a d'ailleurs Djc = D^, on pourra prendre les diffé- 
rences par rapport à z. En faisant, avec M. Tchébychef, 'i, = D*<ï>j, 
le terme le plus élevé de t}>.ç devient 

\)s ^is — {?ll: -.ç _i- 
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Pour réduire le coefficient de z^ à l'unité, nous poserons plutôt 

(54) ''''"' (ï^ D*4>,= C.D^4>,. 
On trouve alors tj/, = A* a?*, et 

ou bien 

, n-4-i , -, . (n-+-i)(n-+-2) 

1^1 = 37 , ^^=z X^ — in-^l)^-^-- '-^ '9 .... 

2 2.3 

I 

Quant à la démonstration de ces formules, elle est facile à l'aide 
de certains théorèmes du calcul aux diflTérences finies. En faisant 

X = D~* on a, comme on sait (* ), 

n n-hl 

/l+/2-+-...-H/«=2//= S/'» 

1 1 

puis, pour deux fonctions/, F, 

[^ Fx • ^/x = ^x -fx — [^ Fa;-4-l • ^/xy 

et par suite 

(55) { '^ 

Toutes les différences Df qui se rencontrent hors du signe ^ se 

déduisent de la série de valeurs yà-, /x^^ , . . . ^fx+s-i ? correspondant 
aux arguments x, x -\- i, . . . ^ x -hs — i, qui, pour les deux limites 
x = ieix = n'\-i, deviennent respectivement i, 2, ..., 5, et 
/i-t-i,Ai-h2, , . ., n-{- s. Or, en vertu de (5 1), ^s s'annule pour 



m 



(*) Dans les formules de M. Tchébychef, 7 signifie aussi somme jusqu'à m 


exclusiçement. 
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ces arguments. Il s'ensuit qu'en faisant /:= ^j, les termes hors du 

signe [^ s'évanouissent aux limites et qu'on aura, d'une manière 

générale, 



ou bien 

n 



(56) 2F.D^*,(a7) = (-i)^2^^(^-^*)^'P» 

1 1 

^s étant défini par la formule (5i). Lorsque F représente une 
fonction rationnelle d'un degré r, inférieur à 5, on a D^F = o, et 
il vient 



rt 



(57) 2f.D^4>, = o. 

1 

Dès lors, en posant F = D''^r> '' <C ^j nous aurons 

(58) ^D''4>r.D^4>, = o, 

et cette relation a lieu, évidemment, pour r^s. C'est précisément 
la relation (3i), ou (48), qui suffit à déterminer les fonctions ^s 
à un facteur numérique près. On en conclut que 

(54) i^,(^)=C.D^4>,. 

Pour transformer les numérateurs (/.^f), nous avons, en vertu 
de (56), en écrivant i pour ^, 



n n 



(59) ^^M. = '^/i.Ds^s(i) = (-iy^^s(.i + s)Ti'fi, 

1 1 

OÙ 

(60) 4>.ç(t -h 5) = i{i -h 1). . .(ï -f- 5 — i){i — n). . .(i— n-h s — i). 
On a 

'î>o = 1, 'l>i(l-{-i) = l{i — /i), *2(^'-+- «)= ^'(^' -\-i){i — n)(i — n-h 1). 
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Il reste à traQsformer les dénominateurs /^4'.?' ^^ faisant, dans 
la formule (59),/= ^j, on trouve 



n 



(61) i2'î''=(-')*2**^'"^*^'^*'''^^'^- 

1 

Or, le terme le plus élevé de ^s étant x% on a D*^Ji^ = 5!, et en 
remarquant encore que ^^(t -f- s) s'annule depuis i 7= n — 5-4-1 
jusqu'à i= n, l'équation (61) se réduit à celle-ci 



n- s 



(62) I^2«^,? = (-iy^!2*-(*'-^^)- 



s\ 



On a d'ailleurs, en faisant toujours C= - — ^? 

(25). 

(63) G.DHi — s). . .(i-H5 — i)= i(i-+-i). ..(i -+-5 — 1); 

en nous servant de cette relation pour transformer <^j, et introdui- 
sant 1^ à la place de \^, (62) peut s'écrire 



n — 54-1 



(64) ! iS^'"^"')'"- S (l-n),,.(i^n-^s-i) 

[ X D*(i — 5). ..(i-1-5 — 1). 

Si nous appliquons la formule générale (55) à l'expression qui 
précède, il vient 



(n—s+l 

( xDs{i — n).,.(i — n-+-s — i), 

car les termes hors du signe disparaissent aux deux limites [ceux 
qui dépendent du premier produit (i — n) ... à la limite supé- 
rieure, à cause du facteur i — /n-5 — i, et ceux qui dépendent 
du second (i — s) ... à la limite inférieure, à cause du facteur 

i—i]. Or 

D*( t — /i ) . . . ( i — /i -+- 5 — I ) = D*' /•* = A- ! 
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et 

tït + i). . .(t-f- 25 — i) = D(i — i)i ... (i-\- is — i). 

V / V / 25-hI V / V / 

En substituant, intégrant et mettant les limites, on trouve 

(66) ^ ^ ^'-^ ,, 

\ {2s-\-i){n — 5 — i)! 

On pourrait établir cette formule à Taide de l'expression générale 

^ 25-M L 2s{n — l) J 

mais la démonstration qui précède peut suffire. En substituant (Sg) 
et (66) dans la formule (46), on obtient la série (52) de M. Tché- 
bjchef, qui devient 



2.3 



/-^^/'-^(^^•«"Sr^'v-H-... 



y;! t i-hs — in — i n — i — s-\-i^ - 

s\n{n^—i)...{n^ — s^)^i si * •" 

à cause du facteur G, introduit par la formule (54). 

27. Il faut enfin chercher une expression de t}»,, appropriée an 
calcul numérique. Nous avions 

(5i) ^s{^) = (a7 — i). . .(a? — s){x — n —i). ..{x — n — «), 

(54) t^,= G.D^<ï>,= ^D^<^,. 

En posant . 

(68) Fs{x) = {x — j)(X'- •!).., (x — s), 

on aurait 

Fs(s^i-x) = i-iyF,(x), 
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puis aussi 

(69) { DFs(^)= -^Fs(^), 

X ^^ S 

^ {x — s){x — s-\-\) 

ou bien 

s' 
D'«F^(a7) = 5(5 — i). . .(« - m -i- \)Vs-m{x) = ^- — — -j F,_,„(a7), 



D'«F,„+,.(;r) = (t + i)...(i-H/n)F,(a:)=^-^i;ï-^F^(:r). 



( t -+- m ) ! 
7!" 



Ensuite 

(ar — n — i)...(a7— Al — s) = F(a7 — /i) = (— i)*F(/H-5-i-i — x) 

et, par conséquent, 

(70) *5(iF) = F,(a7).F,(a7— 71). 

On peut maintenant développer F(^ — n) par la formule des 
différences 

F(j7 — /i) = F(5-M — n)H "~ """ DF(5H-i — /i)-h...; 

d'où, par les relations (69), 

F(a7— n) SX — s — i s{s — i) x — s — \ x — s — a 



F(5H-i — n) I \ — n 1.2 I — n 2 — n 

ou bien 

, . Yix — n) SX — s — I 

; F(/l) I n-l 

(71) l 

\ S(S — l) X — s — IX — s — 2 

( 1.2 n — I n — 2 

et, à cause de (70), 

^^ ) » sis — i) _, , . 

. \ 1.2 (/i — i)(/i — 2) ^ ^ 



• • • 
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Il s'ensuit que 

Or nous avons 

D^F,4-2 = ^i±^ (ar - i){x - 2), 
et , 

(-i)ss\Fs(n) = ^sis-hi); 
par conséquent 

D^<ï>,(a:) = *,(s-f-i) I -— ^; T 

(73) < L »•! 'i — « 

^' J (5 — i)5C5-4- i)(5-f-2) a; — IX — 2 "I 

( 1.2. 1.2 n — in — i"'\' 

ou bien 

5!^! V^-0^^* (5-hX)! (57-0! (n — l — i)\ 



(74)W-)=(^2 



X!X! (5 — X)! (37 — X— i;! (/i — 5 — i;! 



de sorte qu'on peut écrire 



5' ç' 
^s{o^) = (- i)^ ^—p (n - I). . .(Al - s) 

[S(S -{- l) X — I (s l)...(5H-2)a7 — 107 — 2 ,"| 

I ^^ ! ^^ . .. 
I n — i 1 . 2 . 1 . 2 n — in — 2 J 

+ (^-;)---(^ + ^) F,(x)+...1, 

1.2. I.2(/l — l)(Al — 2) J 

OU bien, à rebours, 

i^,(x)=Fs{x)-s^'!^Fs-,{x) 

( f\\ 

J 52(5 — 1)2 n — sn — 5-f-i_, , . 

' H r5_2(a7) — . . . . 

\ 1.2 25 2 5 — I 

On a donc ^/i(^) ^= F„(a:), et i^,i[xi)^=o pour «= 1,2, ... n 

28. En posant 

^s{x)= Âo 4- A, F, -+- AjFg-f- A3F3-1-. . . 
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on a, queJle que soit la fonclion A„ 

Di^,(a7) = A, -+- 2 A2 Fi -t- 3 AsFj -f-. . . , 

D*i];^(a:)= i.î2A2-H-2.3A3Fi-h3.4A4F2-f-..., 
puis 

4;,(a7 -h i) = 4^,(^7) -H D J;^(^) = Bo -+- Bi Fi -h BjFî H-. . . , 
Di^,(a? -f- 1) = 04^,(57)-+- D2 1^5(37) = B, -f- 2B2F1 -h 3 B3F2 -h. .. , 

^,(a:-4-2)=4^,(a7)-H2D4;,(ir) + D2<^^(ar) = Go-f-CiF,-r-G2F2H-..., 

>ù 

Ao-f-Ai, B, = AiH-2A2, ..., Bx= Ax + (X-M)Ax+i, 
BoH-Bi, C, = Bi-f-2B2, .... Gx=Bx-v-(X-t-i)Bx-Ki- 



ou 

Bo 

Go 



On a aussi, comme il est facile de le vérifier, 

\ (a:-i)i|>,(^) = BoF, + B,F2-f-B2F3H-..., 

(70) { 

'' ^ J (ar — 2;(3r — i)^5(a^)=GoF2-f-G,F3-f-G2F4 -+-... 

et, par analogie, 

(a? — i)Di];^(a:) = BiFi 4- 2B2F2 -f- 3B3F3 -h. . . 

(X— l)D^^s{^)= I.2B2F|-+-2.3B3F2 

(ar — i)(a7-2)Dî4;^(a7)=i.2GjF2-t-2.3G3F3 



• • • » 



Ces relations permettent de constater, en premier lien, que 
^s{^) satisfait à l'équation aux différences finies (analogue à l'équa- 
tion différentielle des fonctions sphériques) 

'^^ ( -+- [257-+- I— Al — 5(S -m)] 0^1^5(37) — s(s-\-l)^s(^) = Oy 

qui peut aussi s'écrire 

[(.^ -h i)(^ 4- 2)-- -^^^^1 D24; 4- 2(^ -t- i)D4; = *(* + 1)(4^ -H D^/^ 

Cela résulte de la relation 

1(1 -+- i)Gx- (n - i)(X -h i)(X -+- 2)Ax+2 

•+ (3 — Ai)(X -f- i)2Bx-j-i - s{s -f- i)Bx = o. 

qui se déduit des suivantes 

(X -+- I)2AX-H, = ~-^ ; '- AX, 

n^,\n. (X-M ) (/i — 1) — ^(5H-i) ^ 
(A + i)Bx = —— X ; Ax. 

il A — 1 
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29. On s'assure aussi facilement, en réunissant les termes de 
même rang, que (x — i)^s(^) est égal à la somme 

Il en résulte l'identité 

(78) 4/,+t4-(^-^^ ^ps-^-^ 47rz7+^-» = ^' 

qui peut s'écrire 

^2 fil 5Î 

en faisant toujours z = x • Elle permet de déduire J^^^i 

de i^s et ^s-\ • On a d'ailleurs 

En posant 2-3 = -3|, 2*.<}', = -3,, pour éviter l'introduction des 
facteurs 1^ dans les dénominateurs, l'équation (-jg) devient 

^j+i = zi .^5 — 52 — -T — - z^-i, 

45^* — I 

et nous avons l'algorithme 

4(/i2-4) 

-53 = -Sl-^î 5-T 'Si, 

r -^ r T 9(^'-"9) ^2 

0.7 

On trouve ainsi 

7l« — I 



fc„._ 3.--.3 ^.^ 3(„.-,)( „.-,) _ 

^ /i2 — 7 i5/i* — 23o/i*-»-4o7 

^ 18 1008 

, ^3/i2 — 3i , 5n* — 110/12 + 39.9 , . (/i* — i)(/i2— q)(n« — 25) 

^ 44 17^ '47^4 

3/12 — 43 , i5/i'»— 45o/i2-4-2o5i , 35/1^ — 1645 n* H- 17207/12 — 2720 

^^ ' 52 ' 2288 27456 



- 59 - 

La forme générale de ^s{^) est assez compliquée, en voici les 
premiers termes 

*-='-i^,[--'-N^]'- 

s(s — i)(s — i)(s — 3) r . 5(5— i)-f-3 , 

1.2.2^(25 — 0(25 — 3) L 3 

5s* — 2258-4-375' — 205-t-2in 

+ ^^ J^-*+... 

30. Le terme général de la formule (67) peut s'écrire 

En faisant 

(— i)*4>5(a7-h 5) = 37(374-1)... (a? H-5 — i)(n — 37 — 5-4-i)...(/i — a?) 

__ (s -\- X — i)! {n — x)\ 
~~ {x — i)! (n — 37 — 5)!' 

il devient 



R 



(*') '''^(^^ s\in + sy. Z s\{x-i)l s\(n-a>-sy. ^'^^'^^^ 



1 
i 



OÙ les deux facteurs de D^^sont des coefficients du développement 
de (a + 6)^^"* et de (a -4- 6)"-^. 

En arrêtant la série (62) ou (67) au terme qui répond à 5 = / — i, 
on trouve pour /(oc) une valeur approchée /^'^(^), qui, pour 
x=ij sera/^'^(i). Pour la calculer, il faut faire x = 1 dans^|;,(a?). 
La formule (74)? qui se réduit alors à son premier terme, donne 

(8.) M.)=C.*.(, + .) = (-.>'(^(^^^=^, 

et le terme général de/(i), que nous désignerons par #j, devient 

1 
expression qui peut se mettre sous la forme 

1 
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où A, B, C sont des coefficients du développement de certaines 
puissances de {a -f- b). On a ensuite 

Les différences successives D/(i), D^/(i), . , . s'obtiennent aussi 
facilement. Nous avons d'abord 

et la série (72) montre que le premier terme effectif de D"*^, pro- 
vient de lÙ^'^^Fm^{^)t les termes précédents étant nuls; ce terme 
est indépendant de x^ les suivants ont le facteur (^ — 1) et s'an- 
nulent pour x=^i. Il ne reste donc qu'un seul terme, et l'on 
trouve, pour x = \ j 

ou bien 

) ^'^ ^ m\ (25)! {n — s-i)\ (s — m)\ 

En tenant compte de (81), il vient, pour le terme général de 

(85) D-^f. = ^ 1"-- "' - ')' ^i±:^ §s. 

m\ {n — 1)! {s — m)! 

La somme devra être prise depuis s =^ m jusqu'à s = l — i, car, 
pour s <^m^ T^^^^s s'annule, 'i^s étant du degré s. En désignant le 
coefficient de Ss par (m, 5), on aura donc 

(86) D"'/(/J(i) = (m, m)^,n-him, m 4- i) J,„h-i + . . .-H(m, / — OcF/_i. 

Pour le calcul des coefficients, il est bon de noter que 
(o, s)=.i, et 

(m, s) (s -h m)(s — m-h i) s(s -hi)^ m(m — 1) 

{m — if s) fn{n — m) ni{n — ni) 

M. Tchéb^chef, dans son Mémoire de 1876, donne un exemple 
numérique de l'application de ces formules, qui sont particulière- 
ment commodes. 
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31. On peut enfin, avec M. Tchébychef, chercher l'expression 
de la somme des carrés des erreurs ou résidus A^y. 
En parlant de la formule 

(46) /(^)=^„(;^)î^+,^,(^)?^+... 

et nous rappelant qu'en vertu de (48) nous avons 

24^,.<}'5 = o pour r^s, 

nous obtiendrons S/- en ne conservant dans f^ que les carrés 
des termes successifs de (46), 






ou bien 



(88) ^f.^\lB^^V--f^^ 



24-3 2-^1 

En allant jusqu^au dernier terme 4'//-<î ces formules reproduisent 
les valeurs données de/} et de S/^. En nous arrêtant au terme <}l/_^, 
nous trouvons la valeur approchée y^'\ et la différence de (^lauchy 

représente l'erreur de la formule abrégée. On aura d'ailleurs, tou- 
jours à cause de (48), 

(89) 2(A//)2 _ -^^ + 2^.;^^ ^' • • 

et, en comparant (88) et (89), 

(00) .(../),../,_ [^_..._M-iP. 

On en conclut que 

(91) 2(A/+i/)î = 2(A//)î- t^ii' 

OU bien 

(97.) { _ 9J-4-1 [S/(t-i-i)...(g'-^/— r)U — n...(/i— t'—Z-t-OPT"]» 
L\i\ /l(/l2 — i)...(/i«— /i) 
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On a d'ailleurs, pour / = i , 



et 



A/ = /-^V 



2(A/)« = S/*-^(S/)»; 



maïs le calcul direct de A/ est évidemment plus simpl^e. 

32. Je tâcherai maintenant de donner, en quelques mots, une 
idée des recherches de M. Gram, en conservant toutefois la nota- 
tion employée jusqu'ici. Il s'agit de déterminer par la méthode 
des moindres carrés, à l'aide des valeurs données /<, /21 • • • , fnt 
les coefficients A, B, G, ... de la formule d'interpolation 

/= Aaro-h Bar' -4- Gar"-4-. . . , 

où x^^ x\ x^ ^ . . . sont des fonctions connues de x. Si la précision 
des observations est variable et qu'on désigne par v leur poids 
(l'inverse du carré de l'erreur moyenne), le problème consiste à 
rendre minima la somme 

où y, qu'on peut aussi écrire y*^^^, représente la valeur de y* que 
donnerait la formule, arrêtée à / termes (la somme étendue à 
toutes les observations). Il en résulte / équations de condition de 

la forme 

2p(/-/0^' = o (r=o, I, ...,/-!). 
En posant 

( 93 ) 2 vfx^ = */, 2 (^ a?< xi = pijj 

elles deviennent 

(94 ) si= Xpio -f- Bpii + G/)/2 + . . . 

et elles déterminent les / coefficients A, B, G, . , . dont les valeurs 
dépendent du nombre /, tandis que les quantités Si et pij en sont 
évidemment indépendantes. 

Si nous prenons un terme de plus,y*^ devient /'■•■*, et les nou- 
veaux coefficients sont donnés par l -\- i équations de la même 
forme, dont nous retrancherons les précédentes. En désignant les 
différences des coefficients respectifs par ^, i, c, ... et faisant, 
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pour simplifier l'écriture des formules, /= 2, nous aurons les deux 
systèmes d'équations 



So = A/?oo+ B^oi-+- G/>oî, o = apoo-h bp^x -h G/?o2» 

*i = A/?io-4- B/?i,-4- G/?i2, 0= a/?io+^/?ii-h GjOij, 

*î = Ajt?2oH- B/>2,-H Cj922; /(3)_y(2;= «a;o ^_ 6^?' _|_ c^r'. 

En égalant les deux valeurs de C, fournies par le premier et le 
second système, et qui sont des rapports de déterminants, il vient 



d'où l'on lire 



Q _- [/^OO/^ll^î] __ r/^00/^1ll r ^{3)__ A2)l 
[/>00/?lljD22] {PQQPna^'] ^ '^ ^' 

y (3) _ y(2) _. \P^^P\\S<t\ [/^OO/^lia^] ^ 
[/?00/?ll/?22] [/?00i>ll] 



Pour /= I, on aurait 



•^ />00 



et l'on trouve ainsi, par additions successives. 



/?00 



[/^OO/^ll] /> 



[/?00/>ll/?22| VP^^P\\\ 



On obtient donc pour / une série ordonnée suivant des poly- 
nômes qui sont des fonctions linéaires de x^^ x\ x^\ .... En les 
désignant par ^m(^)> nous aurons 



(96) 



'tmCa?) = [/?ooi?ii . . . a:'«] = 



7>io P\\ 



P\m 



• ■• ••• ••• ••• 



x^ x' 



X^n 



et, par suite. 

Enfin le dénominateur du terme général peut s'exprimer par S(^<{/f,^. 
En effet, on voit facilement que ^v^m^r s'obtient en remplaçant 
^', dans ^mt par ^^oc^^m\ or, cette somme est nulle pour i<Cr^ 
car elle s'obtient en remplaçant, dans ^r? la dernière ligne x^, 
x\ ..., x'^ par /?/o. Pis, . . ., pir] pour « = r, elle représente le 
déterminant 

Il s'ensuit que ^r'i„i'ir s'annule pour m < r, et, par raison de 
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symétrie, aussi pour r <C m\ pour /• r= m, il reste le produit de Pf'"^ 
par son premier déterminant mineur. On aura donc 

\ ^v^,''^s=o pour r^s. 

Ces relations expriment la propriété caractéristique des fonctions i^ 
qui peuvent servir à développery* en séries selon le principe des 
moindres carrés, La « série d'interpolation » a la forme 

OÙ <|/o (•^) = ^^' En l'arrêtant à ifmt elle représente le résultat de la 
compensation des observations données par une série procédant 
suivant les fonctions x^^ x\ . . ., x^. C'est proprement une for- 
mule d'interpolation si le nombre des termes (m + 1) égale n, et 
une formule d'approximation pour m-\- \<^n, 

La série (98) coïncide avec celle de M. ïchébychef dans le cas 
où j:®, x\ ... sont les puissances de x. Lorsque x^^ x\ ... sont 
des fonctions i qui satisfont aux conditions (97), on a évidem- 
ment pijT=o pour i^jf et les coefficients A, B, C, ... prennent 
la forme caractéristique de la série (98); on a, par exemple, 



// H 



/?22 Zvor"x' 

On trouve enfin, comme on le vérifie aisément, pour la somme qui 
devient un minimum, 

(99) M/ = 2t/2_Q,, 

ou 

et M/ tend vers zéro à mesure que / augmente. 

33. M. Tcliébychef n'avait considéré que les séries d'interpo- 
lation qui découlent d'un développement suivant les puissances 
ascendantes de x^ en particulier pour p = i ; dans ce cas, pour des 
ordonnées équidistantes, la série (98) donne nos formules (Sa) 
ou (O7). M. Gram indique encore une autre série pour v =z nx* 
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En prenani di/= - — > on arrive à la formule d'interpolation 

(a? Xi) 

de La grange 

Les séries trigon orné triques par lesquelles on représente une suite 
d'observations d'un phénomène périodique (Bessel, 1827; Astr. 
Nachr., n® 136) fournissent un autre exemple de séries d'interpo- 
lation, construites selon le principe des moindres carrés. 

Lorsque la série doit servir à représenter une fonction continue 
/(^r), on remplace les intervalles de l'argument x par dx^ les 
sommes se changent en intégrales définies, la première équa- 
tion (97) devient 

(100) I V^r^sdx =0, 

et le développement a la propriété de rendre minima l'intégrale 



f v{f-fiydx, 



p OU f^^^ étant toujours la valeur de y* fournie par / termes de la 
série, et v une quantité positive fonction de x. On se trouve ainsi 
conduit aux. fonctions V, étudiées autrefois par Sturm et Liou- 
ville(i836-i837). 

De ce genre, sont les séries de Fourier (*), qui procèdent sui- 
vant les sinus et cosinus de mx\ on a ici ç^ = i, les limites des 
intégrales sont -h u et — t:. 

Les fonctions sphériques, ou polj^nômes de Legendre (2), satis- 
font à la même condition en prenant ç^ = 1 et pour limites +1 
et — I. De même, les fonctions de Bessel, et d'autres polynômes 
qui sont étudiés dans le Traité de Heine. M. Gram, dans sa ïhèse, 
a indiqué une grande variété de séries d'interpolation; il étudie, 
notamment, celles qu'on obtient en partant des formules 

f—{a-^bx-^cx^-\-...)e-^^ s> — e^ 

(*) Voir la Notice historique de 13. Rieniann {Bull, des Se. math., t. V, juil- 
let 1873), et une Note de M. Tripler {Bull, de l'Acad. de ViennCy 1876). 

(») G. Plarr, Comptes rendus, mai 1857. — E. Rouciie, Journ. de l'École 
Polytechnique, i858. 

R. 5 
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liellemenl aux cas où les données sont 1res nombreuses, où, par 
conséquent, la méthode des moindres carrés entraîne des calculs 
d'une longueur disproportionnée avec les avantages qu'elle procure 
sous le rapport de la précision du résultat. La méthode des moin- 
dres carrés resterait réservée aux cas où le nombre des données 
n'est pas très grand, parce qu'elle permet d'épuiser, en quelque 
sorte, les ressources qu'elles offrent pour diminuer l'erreur moyenne 
à craindre, et aussi pour nous rendre un compte exact du degré 
d'approximation qui a été atteint. 

La méthode de M. Tchébychef a fait, dans ces derniers temps, 
l'objet d'études importantes, publiées par M. O. Backlund dans 
le tome VII des Mélanges mathématiques de l'Académie de 
Saint-Pétersbourg (mai 1884) et par M. P. Harzer, dans les 
jjos 2757-2759 des Astronomische Nachrichten (novembre 1886; 
voir aussi BiilL astr,, t. III, p. 61 1). Ces études, qui sont 
accompagnées d'exemples numériques, nous faciliteront notre 
tâche. 

35. Il s'agit donc de déterminer les coefficients A de la formule 
d'interpolation 

(i) /(^)= Ao-H AiiP-H.. .-h A„a7« 

par des relations de la forme 

c,A| = Si/— Si/+ S3/— . . .± S^/, 

où S|, S2, S3, . . . sont des sommes de valeurs successives de/, 
ces valeurs étant données en nombre infini pour des valeurs équi- 
distantes de x. En faisant 

(2) c/ = (Si— Sj-h Sg — . . .)^'î 
on trouve d'abord 

(Si— SjH- S3 — . . .)/= AoCq-i- AiCi-4-. . .-+- A^C/j, 

et, cette expression devant se réduire à un seul terme A/C/ (avec 
un facteur ci aussi grand que possible), on a les conditions 

(3) Co = 0, Ci = 0, ..., C/_i = o, C/^.l = 0, ..., c« = o. 

Si les valeurs données sont très rapprochées entre elles, les sommes 
pourront être remplacées par des intégrales ; nous pourrons encore, 
sans nuire à la généralité des formules, supposer l'argument x 
contenu enlre les limites -h i et — i. Cet intervalle étant divisé 
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en une série de Ironrons, ayant pour limites successives 

— I, rtj, «2» •••) <^/«» +ï> 

et correspondant aux sommes partielles S|, S2, . . . , si nous repré- 
sentons par le symbole / la suite d'intégrales définies 






ou bien 

pftl jr^ftt >,«3 ^rtj /*"*"* 

«-^—1 «^rtj *^a. «^«4 »>'«,„ 

nous aurons 

%J m 

F(^) étant l'intégrale indéfinie f /(x)dx. Ensuite 

(4) CiAi= I f.x, 

OÙ le second membre s'obtient par voie de quadrature mécanique, 
à l'aide des valeurs données de/, tandis que le coefficient c/ ou cj"*' 
peut être calculé, une fois pour toutes, par la formule 



m 



(5) 


C^m) __ 


j x^ dx^ 


d'où l'on tire 


jn*^ 




(f)) '■^^c^'^J-a'i+i 


a^f^-\-.. 





Les limites intercalaires a^, «2» ...,«;» se déterminent par les 
conditions (3), en égalant à zéro l'expression (6), dont le dernier 
terme est, tour à tour, o ou zh i. 

En prenant m = /i, on aurait donc n équations (non linéaires) 
pour déterminer n inconnues (a<, aa, ...) pour chaque valeur 
particulière de i. Il faudrait ensuite, parmi les racines réelles, 
choisir celles qui forment une suite régulière entre les limites — 1 
et -h I, et qui donnent pour Ci la plus grande valeur possible (*). 

(») On ne peut guère recourir ici aux équations de condition û?c„= o, rfc, = 0, ... 
que fournit la théorie des maxima relatifs; elles seraient en général incompatibles 
avec les équations données, à moins de prendre m > /i et d'admettre des racines 
égales (n» 40). 
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Le mode de solution employé par M. Tchébychef, fort ingénieux, 
mais d'une applicalion assez pénible, repose sur les remarques sui- 
vantes. 

On a, en développant suivant les puissances descendantes de x^ 

[ I ( — 1)« 2 



y x-\-\ X — I X — ai X — a^ 

à cause des relations (3), en vertu desquelles les coefficients ci 
s'annulent jusqu'à l^=n^ excepté pour l=.i. En multipliant 
par dx et intégrant, jusqu'à ^ r= oo, on trouve 

, a? -f- 1 (x — «2)2... Ci Cn-h\ 



Cette relation peut s'écrire 

fo. ^ + ' (x-a,)Hx-a,y,.. ^ -£ù^££±i^£n±.v... 

^^ (a: — 1)^1 (ar -ai)2(a7 — as)*... 

On pourrait aussi la déduire de la formule générale qui exprime 
la dépendance mutuelle des coefficients d'une équation algébrique 
et des sommes des puissances semblables de ses racines : 

•^i Si .ï, 



(9) ix — a){x — b){x — c).,.=x^.e -^ «x» s.» '-'^ 

36. M. Tchébychef se sert de la formule (8) pour déterminer Ci 
et les racines a^ a^, . . . , à l'aide de raisonnements empruntés à 
la théorie des fractions continues, en considérant les réduites qui 
se rencontrent dans le développement de l'expression 



/ 



(x — i}^i 

X -^ l 



eir 



«+I 



Mais il convient peut-être de circonscrire tout de suite le pro- 
blème, en constatant que la résolution directe des équations (3), 
en supposant m = ai, donne toujours n racines symétriques (c'est- 
à-dire égales et de signes contraires) dans le cas de i=L n^ n — 2, 
n — 4? • • • et Al — I racines symétriques (la /i^'""*^ étant ± i) dans 
le cas de f = Ai — i , ai — 3, . . • . On peut, en eflct, prouver que la 
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soliilion est fournie par ai ou ai — i racines symétriques dz a, di 6, 
zhcy ... (dont Tune égale à zéro, si leur nombre est impair). 

Pour le démontrer, remarquons que, si nous changeons x 
en — Xf la formule (8) devient 

x-^i (x^a,)Hx-^a,y.., ^ ±-£;-±^;-H... 

Lorsque i est un nombre pair, en multipliant (8) par (lo), c/ s'éli- 
mine, et l'on trouve, suivant que n est pair ou impair, 

i' . x^ — al x^^al %-! 

\^P^'' :;;¥— :;lzrr^---=^"""' 

(il) i pair 



x^ — a| 


x^ — al 


x^ — af 


x^ — a| 


X^ — I 


x^ — al 



fn-^i ^ 



n impair. . . —g r- — r ^ . . .== er'»^' 

^ a?* — af a?*— a| 

Lorsque i est impair, on élimine Ci en divisant (8) par (lo); on 
trouve ainsi 



(f2) i impair 



X -h ai X -h a^ X — a% X — a^ —h-*- 

n impair. . . • • • . . . = e^ ^ 

X — a\ X — as /*. -1- ^- /*• _L_ -T». 



x-\-\x-\-axX-^a% x — a^x — a^ ~xî+ 

n pair • • • • . .= &^'"'^* 

^ X — \ X — a\ X — as a^-hajiCH-a^ 



X 




a% X 


— a^ 


X 


-f- 


a^ X -\- ai^ 


X 




a\ X • 


— a,, 



Comme le degré des dénominateurs lie dépasse pas n ou /i -h i , 
nous pouvons remplacer les exponentielles par l'unité. Les rela- 
tions (il) signifient dès lors que les quantités aj, «3, ... coïnci- 
dent avec «2» «^, ... ; dans le cas de n impair, il y en a une qui 
est égale à l'unité. 

Or, on ne peut faire a< = «2? . . . , car la formule (8) conduirait 
alors à des conditions incompatibles; il s'ensuit que les racines ai, 
«8, . . . sont égales à — «2? — ^4> . . . , en d'autres termes, que les 
racines sont symétriques (dans le cas d'un nombre impair, il y en 
a une qui se réduit à l'unité, de sorte que m = n — i). 

Les formules (12), relatives au cas de i impair, nous disent que 
les racines sont encore symétriques, mais l'une se réduit à zéro; 
et, si n est pair, une autre se réduit à l'unité. Pour le comprendre, 
il suffira d'écrire les équations pour ai = 4 et ai = 5 

/! = ♦ 

{x -f- i){x ■+■ ai)(x -h a^){x — a^{x — ai^) = {x — i)(.r — - ai)(a? — «sX^ -H <^2)(^ "*" ^4)» 
(x -f- aiXx -+- a^Xx H- a5){x — a^X^ — a;) =z{x — ai){x — a3){x — as) (a: -H ai){x -4- a^) 
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En faisant a:=:o, elles donnent respectivement r7|^2^3<^i = o et 
ai «2<^3<^4Û55 = 0, d'où il suit qu'il y a toujours une racine nulle. 
En faisant, dans la première, x = i^ le second membre s'annule, 
et il s'ensuit qu'il faut faire ^4 = 1 ou «i = — i , ce qui réduit le 
nombre des racines utiles à trois. On a donc encore m = n — i . 

37. En résumé, le nombre m des racines est égal à /i, lorsque n 
et i sont de même parité (i=zn^ n — 2, n — 4> •••)» ^° ^"^^ 
m = n — I (l'une des racines coïncidant avec ± i)pour i=z n — i , 
n — 3, .... Nous n'avons plus, dès lors, que des racines symé- 
triques dz a, d= 6, ±Cy . • . (y compris la racine zéro, si leur 
nombre est impair). La formule (6) devient 



03) \±ic.^m.^(-^y-^(-^)- 



(J-«'-^»+6'^*-...), 



et Ton voit que, de cetle façon, on aura, pour m = n^ 

C,i—i ^ O, Cfi—3 = 0, . , . , Cn-k-i = O, C/14-3 = 0, . . . , 

et pour m = n — i, 

Il reste donc, dans le premier cas, les conditions C/i = o, 
c,i_2 = o, . . . , et dans le second, les conditions c,,_« = 0, c„_3 = o, 
qui, les unes et les autres, peuvent s'écrire 

(i4) \ 1 = a'»-ï— 6'«-i-hc'«-i — ..., 



Le nombre des inconnues a, 6, c,'. . . est — > ou bien — ; — (plus 
la racine o). En même temps, 

i = m, m — 2, m — 4> ••• cl n = m ou m-hi. 

11 est entendu que celle des équations (i4) qui correspond à c/=: o 
doit être remplacée par la suivante 



(l5) - — (~I)' ^-^ Q:r3a' + Ï— ^>'>i+C'»'— .... 

2 4 



— 72 — 
On trouve ainsi, par exemple, 

m = 'i (i = *i ... «=}, \ — a^=lCi, 

'^ = 3 l t = 3 ... aî = i, a* — t = C3, 
/i = 3 ou 4 { 

i f = I /iV :^ 1 /ï2 1 — 1/1. 

6=0 I*—» ••• «—2» « J — ï^l» 

de sorte que les limites intercalaires et les coefficienls c/ deviennent, 
pour w = 3, 

Limites. Coefficients c. 

1 



,x 3 



Co = — 1.1748 



2^ 

.1. 

^^3;±(^) et o c,=: 4-0,4142 



Cz — -HO, 5 



2 

X 
— [-) et o C3 = — 0,25. 



On a ensuite, par exemple, pour £ = 2, 

_! 1 

7,^i=( fdx— I fdx-\- f fdx, 

2 2 

oïl les intégrales devront être évaluées par quadratures méca- 
niques. Pour m = 4? 5, 6, . . . , on trouve de même 

1=4 ... « — ^ = à^ — 63= -, à^ — ^5= c 

1 24*' 

'^~"* /4rr=2...a— ^> =«5—65= 1, «3 — ^>3 = 1 _ 5 c„ 

/i = 4 OU 5 \ îi 24 



i = o ... «3 _ ^3 — ^5 — b^— -y a —b = Cn, 

2 2 4 "' 

t = 5 ... a2— ^>2=a^— ^>i= i, a6 — 66=1±^, 
- . 2 2 

m = 5 

AI = 50U6 / t = 3 ... a2--^>î=rt6_^6= 1, «4_Ô4= 1±_^ 

] • 2 ' 

c = o J 

t=i ... a*— 6*=a6— ^>6= 1, a2-62=li-^, 

•2 2 

1 = 6 ... a — 6 H- c = a3 — ^>3 _^ c3 = a5 — ^5 _+. cs = I , 
m = 6 ) ^' 

/i = 6 OU 7 i ^' ~ ^' "^ "'^ = 2 "" 4 "'fi' 



- 73 - 

38. La résolution de ces équations n'offre aucune difficulté 
sérieuse. 

Dans le cas particulier de i^ m, ainsi que l'a remarqué M. Tché- 
bychef, les valeurs des constantes peuvent s'obtenir d'une manière 
explicite. Les conditions sont ici 

1 4 

(•6) ^ . 1 

' a'"-»— 6'«-*H- c'«-»— ...= -, 

2 



« 



On y satisfait en posant 

a = cos » «> = cos > c = cos 

m -4- 1 wn- 1 m -h I 

Cela résulte des relations 

iir 2*^ , /iir î 

cos cos h. . .dt cos = - , 
2AI-+-I 2/H-l 2AH-I 2 

J 11 21: , /i — T I I 

I cos— — cos — -t-... ni cos ir = - in -, 

\ n n Al 2 2 

comme on le vérifie aisément. 

Les formules (3) et (i 3) montrent que, pour 1 = m, nous avons 

On trouve ensuite 

(20) c'.o'~ — ^— - ^(1) - """ 



*^* — T7r~T~T ' *'2/*-Hi — 



Q/IH- I 2/1-1-2 

et, pour m ;> I, 

i ^r/.) _ (-iV^» m-f-i mH-5 mH-6 

ou bien, pour m >2/, 

(22) ( — I )'« 2'«"» ci;«-2''> = mi — /7Î/-1, 
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où le second membre représente les coefficients du développe- 
ment de (i — a:)(i + x)"^. 

39. Dans le cas général, on peut se servir avec avantage des 
relations bien connues qui existent entre les coefficients/?,, d'une 
équation algébrique et les sommes 5,, des puissances semblables de 
ses racines, à savoir 

ou bien, si l'on fait s" — Sn = nçny 

Ç3 = PiPi—P3i q& — {p]^Pi)Ç3-^PiP!t — p5^ ' 

^7 = ip] — pi)qs-^ {piP3—p%)qz-^pip&—pii — 

On trouve ainsi, par exemple, 



m. t. 



I 



4 4 . . . 4 « — I = 4 ^ H- 1 = /5 ; C4 = - • 

o 

4 2... (4a — i)s = (46-f-i)2=2/Ï3 — i; 4^2=3 — /^. 

5 

4 o... i{a — by — 5(a — ^)3-h 9(/z — b)=-< 

5 5 . . . a = - v/3 , b = -; Cs = -' 

2 2 16 

5 3... 4«' — I = 46*-f-i = /S; 4^3= /5— 2. 
( (a« — ^»2)* — 2(a2 — 62)-H 7=0; 

^8— «6— -a*-h-=o; ) =y-+-ci. 

3 2 4 \ 2 / 4 

loi I 

a3 «2 a + o = o ; 

228 

TZ j ITZ StU I 

a = cos-> 6 = cos — > c = cos — ; 0^= ir- 
1 7 7 32 

I a3~ i aï— g (/^ — i)a + -jîg (/I5 - 3) = o; 
( C4 = — -j^ (/iÏ9 — 5). 



5 I. . . 



6 6. . . 



6 4 



I a'— - a*+/)îa— />3 = o, 

• • • ■ 



/'i + Jy'l-J/'.-^' </'»+87;; + '="- 



75 — 



m. 



6 






3i5 



Pî 



175 



jOjH- 18 =0, 



Ç& = (/^î '—P%)ÇZi qz=P\Pi —Pi 



["'-'v 



«2 H- C* = I , 



a* — a* = c' — c^=-> û?=o, 

o 



a = cos-> 

o 



4 



sin 



TT 

8 



^' = -64' 



Voici d'ailleurs les valeurs numériques des constantes 

Limites intercalaires. 



m. 
O 

I 

2 

3 

4 



n. 
OOU I . 



Coefficients. 
Co=a 



4ou 5. . 



I ou 2. . . Ci = — I 

20U 3. . . Cj = -f-o,5 

Co = — 1,1748 

3ou4... C3 = — 0,25 

Cl = -f-0,4l42I 

C4 = -f-0,l25 

Cj = — , 1 5 1 39 
Co = -+-0,83443 

C6 = — 0,0625 
C3 = -1-0,05901 

Cl =—0,22777 

60U7. .. C6 = -}-0,o3l25 
C4 = 0,02407 

C2 = -ho,o6653 
Co = — 0,64758 

7 ou 8... C7 = — 0,01 562 



o 



5 5ou 6. . 



d=o,5 
±0,79370 

±0,70711 eto 

±0,84090 eto 

± o , 80902 ± o , 30902 

±0,87305 ±0,37305 

±0,89729 ±0,60591 

±o,866o3 ±0,5 

=t 0,89945 ±0,55589 eto 

±0,91681 ±0,67411 eto 

±0,90097 ±0,62349 ±0, 22252 

±0,91975 ±0,66389 — OîM4i4 

±0,93091 ±0,72074 ±0,29483 

±0,93839 ±0, 76036 ±0,48387 

±0,92388 ±0,70711 ±0,38268 eto 



et o 



Pour mieux faire saisir l'esprit de la méthode, transcrivons ici 
les formules finales pour /i = 3. 

— i,i748Ao= / fdx— j fdx-^- I fdx\ «=0,7937, 

o,4i42Ai= / fdx — I fdx -h j fdx — / fdx; « = 0,8409, 
jAj= / fdx— j fdx -h j fdx; a = o,5, 

X» — rt >^0 x,rt ^1 

— {^3= 1 fdx — / fdx-^- I fdx— I fdx; « = 0,7071, 
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Les intégrales s'obtiennent par quadratures mécaniques, à l'aide 
des valeurs données de y, que l'on suppose nombreuses et rappro- 
chées entre elles. Il faut, bien entendu, commencer par ramener 
à dz I les limites Xq, x,i de l'argument donné x, en posant 

/_o\ Xw -f- Xo , X„ — \q 

(23) X = ! X, 

•2 2, 

40. On pourrait aussi établir des formules analogues en prenant 
m^n. Ainsi, pour /i = 2, /n = 4? on aurait les limites di a, ±: 6, 
avec les conditions 

la meilleure valeur de c^ s'obtiendrait pour « = -f- j, 6 = — |, 
à savoir Co == |. On trouverait alors 

44 44 44 

Mais cela revient à introduire encore des facteurs numériques 
qui multiplient les intégrales, et ce serait déjà une extension du 
problème. Il est vrai que la relation 

|A2= r fdx — ^ r*fdx-\' Ç fdx 

~4 4 

ne serait guère moins simple que celle-ci 

1 1 

^As = / fdx — 1 fdx-\- j fdx. 



41. M. Tchébychef a fait l'application de ses formules à un 
exemple numérique, emprunté aux recherches de Kopp sur la 
densité de l'eau. M. Backlund a traité, à son tour, quelques 
exemples de nature à intéresser les astronomes (réduction d'ob- 
servations spectrales, comparaison de deux Catalogues, distances 
mesurées d'un couple). 

L'évaluation des intégrales qui entrent dans ces formules peut 



— n — 



élre abrégée par remploi d*un procédé graphique; on peut aussi 
commencer par régulariser les données d^observalion en formant 
des « lieux normaux » de même poids. M. Tchébvchef a d^ailleurs 
donné une formule très simple pour le calcul approximatif des 
intégrales, qui se déduit de la formule des trapèzes, et que Ton 
peut, avec M. Backlund, ramener à cette dernière par Tintroduc- 
lion des valeurs de f{x) qui correspondent aux limites interca- 
laires. Supposons que la limite a tombe entre les valeurs Xj, x^ 
de l'argument pour lesquelles on connaisse y^^ /"<, on fera sim- 
plement 

Xj^ — Xi 

On aura ensuite, par la formule des trapèzes, en supposant 
^4,^:5, . . . , X9 compris entre a et 6, 

_ 5^4 — a ^ Xi — a ^ Ta — .r^ 3^9 — a?? 

— Ja -' /* H /* ■+" — "J /« 

ÎA 2 2 2 

2 . 2 "^ 

(xi — a)*/»— ^^a — «)*A . Xi—^t ^ , ^« — X* ^ . 

2 •" 2(3-10 — 3:9) 

42. Un des inconvénients de cette méthode, c'est qu'il faut 
fixer d'avance le degré de la formule d'interpolation, et recom- 
mencer tous les calculs, si le résultat n'est pas satisfaisant. M. Tché- 
bychef a cherché à remédier à ce défaut par sa seconde méthode, 
fondée sur la formule 

(M) /(^) = oo^o-*- ajOiH-. . .H- a«6rt, 

où 9o> ^M • • • sont des polynômes donnés de la forme 

(25) 0/= ^x^-h^xf-^-h..., 

tandis que les coefficients ao, ai, ... se déduisent des valeurs 
données de /en faisant 



(26) a 



= f/{x)dx, 

*J En. 



m 

m 
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les limites intercalaires étant ici exclusivement celles fournies par 
les formules (i6), de sorte que 

«/Q •^— l •'1 *^— 1 «-0 *^J *^—\ •^_* ^A-^ 

2 i 



• • ■ • 



D'après ce qui a été dit, ces limites vérifient les conditions 

V-*// ''0 — "> ^1 — ^' •••> *'//i-l — "» *'//t+l — "> ^/«+3 — "> •••1 
OÙ 

On peut déterminer les Q de plusieurs manières. Remarquons 
d'abord qu'en substituant l'expression (24) dans (26), on doit 
trouver une identité, d'où il suit que 

/ 0,;, rfiF = I , 

, ^ m 
(28) 

f l ^idx =0 pour / = o, I , . . . , m — i , m -h i . . . w. 

Si 9/ a la forme (25), ces conditions sont déjà remplies pour / = o, 
I, . . ., m. — I, m -4- I, m + 3, . . ., à cause de (2^); il ne reste 
donc que les conditions relatives à /=m, m +2, m + 4, •••, 
qui suffisent pour déterminer les coefficients de 6/, puisque 






On aurait, en faisant successivement m = l^ l — 2, / — 4? • • • > 

I = Pc/', 

les coefficients c^"'' étant donnés par les formules (19). On trouve 
ainsi 

Oo=-j, Oi=— a:, 02=7.37*— 5, 0»= — 4a73-h9.iP, 
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11 est d'ailleurs visible que la forme (25), adoptée pour 8/, est 
exigée par les conditions (28), qui montrent que 8/ est nécessai- 
rement une fonction paire ou impaire, de degré /. 

43. Nous allons maintenant déterminer les polynômes d*une 
manière plus directe. Il convient, pour cela, de poser 

(3o) f f{x)dx = ¥{x) 

d'où 

a„i = [F]=-F(-i)-t-9.F(a,)-2F(a2)-4-...-f-(-i)'«F(i). 



tn 



Nous ferons en même temps 






(3i) i{m-\-i) I 0,„c/a? = (— i)'"X,„+,, 

de sorte que (24) devient 

{2F(a7) = aoXj — ^aiX2-h^a2X3 — . .., 
2F(^)=X,[F]-iX2[F]-hiX3[F]-.... 
d 1 2 

En désignant par z une quantité arbitraire, et prenant 

f{x)= — ^-— , F(a7)='*^r— î^i, 

on aurait donc 



(33) 2^ 





\_z — x\ yz — x\ 1 \^z — x\ 



Or, d'après ce que nous avons vu, par les formules (6), (7), on a 
(pour iz= m) 

\_z — xj z-hi z — ai 2 — «j s — I 

(m + i)ci,?' , (m + 3)cai . 

H :r~r~. -t-. . . . 



m 



J//H-Î 5"*"*"* 



Il s'ensuit que {m -f- i)c5^' représente le numérateur de la fraction 

> et qu'on aura 

z — xj ^ 



[z-x\ (z^-i 



(m -hi)c 



{m) 



^^'^^ '- -' ' "' ■){z-a,){z-a^),,.(z-a,n) 
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Or 



c;;;" = (— 1)"»2»-'«, 



7C aTT 



et les racines a», ao, . . . sont cos > cos ? • • • • En faisant 

5 = cosJ^, et nous rappelant que 

sin/iÇ = 2«-isinu cosÇ — cos— j( cosÇ — cos-^ )•••( cosÇ — cos ir j, 

il vient 

r I 1 (-iV/«-^i.9.(m-4-!) 
^ ^ „jL>5 — x\ sinÇ.sin(m -h i)Ç 

Si l'on fait encore 2>3 = ^ ~h tj» on trouve 
et, par suite, 

En substituant cette expression dans (.^3), et mettant m pour 
m 4- 1 , on trouve 




ou bien 

Si Ton fait x = cos^^, cette relation devient 

(38) 2 T^' = 7^ + ^^"' ''°'"y'' 

d'où l'on lire, en différentiant, 

(39) 2(- 0- ^P"!r^ = - Smp- ^iîV^, 

V î^/ V / I — pa/'* ^ .281117 



puisque 



I /fX,„4_i 



/// -t- i 2 <^/a7 



-Bi- 
ll en résulte les valeurs suivantes des X : 

Xi = COSJ^, X2 = C0S2^ -4-1, 

X3 =cos3j^ — cos^, X4 = cos4j^ — I, 

X5 = cos5^ — cosj^, Xe = cos6j^ — cos2^, 

X7 = cos7^ — cos^, Xg = cosSy — i, 

Xg. = cos9^ — cos3^, Xio = cosio^ — cosa^, 

Xii = cos I ly — cos^, X12 = cos 12^ — cos4j^, 

Xi3 = cos l3jK — COS^, Xi4 = cos l^y — C0S2J^, 

Xi5 = cosi5j^ — cos5jK — cos3^ ■+• cos/, Xie = cosi6y — i. 



On aurait également 



. ^ sin 2 V 

20o=I, 261 = ^^^, 

smy 

-. sin3v i f. sin4r 

sin^ 3 s\ny 

^ sin5r i ^ sinCv i sin2y 

264= . -^ — -z, 205 = r-^ -h - . -^ , 

smy 5 siïi^ 3 smy 



ou bien 

60=1, 61 = — cos^, G2=cos2jK + 3, 03 = — cos 3^ — cosj^, .... 

Il est donc facile d'exprimer les X et les 9 en fonctions de ^ = cos y 
ou de Ç = 2X. 

Ai. On pourrait aussi procéder comme il suit. Supposons que 

Ton ait /(^) = Aq + A^ .r + . . . ; en substituant dans les coeffi- 
cients oLmy on aurait 

Ao-i- Aia^-h. . .= 20/ / (Ao-h Aia7-f-. . .), 

par conséquent, à cause des relations (27), 
(40) ^'« = 0,„c^^'' + 6,„_2ci;r2^-i-. . . . 

Avec les valeurs (21) des coefficients c^'"^, on trouve 

il' , ^ m — I- m (m — 3),, 

1 ' I 1.2 

/7^(m-I)(y/^— 5) '^Z^' . 

1.2.3 //l 4- I 

R. 6 
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où il faut écrire 8, au lieu de Oq, si m est impair. L'intégration 
donne ensuite, en mettant n pour m + i et Ç pour ix^ 



(42) ll^n=,Xn-\-n\n-'i-^— ^^ X;,_4 + . . .+ 2«~i X, (ou a^-sXj ). 

2 1.2 



En remplaçant X| par x ou X2 par 2;r2, on trouverait, pour n 
impair ou pair, 



(43) i ='"-;;""-"M=x„+«x„_,-H^i(:^)x„_.+.... 



Ces formules conduisent aux relations suivantes : 



H 


= x, 








\v 


-X, 


+ 4X, 






H' 


-Xs 


-4-5Xa+i6X, 






ir 


= x, 


+ 7X5+21X3 


+ 


64 X, 


1 
2 


= 8oi 


i$' = e,+ 


2 
3J 




i^ 


= 6». 


H- 36, + !, 







i$« = e«+504 + 99, 



3 2 

7 J 



i?' = x, 

J^ = Xv + 4X, 

||8 = X« + 6X4+i6Xj 

f$» = Xs + 8X6 + 28X4 + 64X1 

^= = 65+463-1^ 

H' = 07+6e5 + l493-8$, 



d'où l'on tire encore 



2X4 = I, 

aX, =$«, 

aX, =? ($'-4), 

2X4 =|»(5»-4), 

2X, =? ($î_4)(Çî_,), 

2X, =$H$'-4)(?>-2), 

2X, =?($'-4)($'-i)(5'-2), 

aX, =^«(^i_4)(5î-2)S 

2X, =$ (5»_4)(5»_3)($«-i)», 

2X,« = $«(l*-4)(?'-3)($»-2)($«-i), 

2Xn = î ($*-4)($»-3)(?«-i)(^-3$« 



•). 
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En développant, on aurait 



•2X1 


-?, 














2X2 


-V, 














2X3 


-|3 _ 


4$, 










t 


2X4 


= 5* - 


41*, 












•2X5 


-?» - 


5Ç5 


+ 4$, 










-«Xe 


-|« - 


6|* 


-+- 8?2, 










•2X7 


-r - 


7?' 


-m4Ç'- 


8$, 








•2X8 


_ç. _ 


8?» 


-4- 20$* — 


i6$« 


> 






2X9 


-$' - 


9$' 


-4-27$5__ 


3l^3 


-4- 


I2Ç, 




'2X10 


= $.«_ 


I0$« 


-h35Ç6 


5o^ 


-+- 


24 S^ 




^Xh 


= «"- 


•II$9 


+ 44^- 


77?^ 


H- 


55^3 _ 


-12$, 


2X12 


-1"- 


I2$l< 


>-4-54$8- 


112^® 


-+-104$*- 


-32e^ 


2X13 


_$n_ 


i3|n 


+ 6j$9- 


i56?^ 


+ 182^5- 


-91?3-M2$, 


tra 1a 


c vqIaiii 


ne Ha 


ft_ ATI ri 


ifTi^rA 


nti 


anl At 





/?l 



260 =1, 
î^e, = $2 -1, 

-263 =$3 -2$, 

264 =:^ - 3$2 +1, 

-265 =$« - 4?^ -t-lï, 

— 2O7 =$"^ — 6$s -f-iof3_ 4$, 
20» =$« - 7$^ +i5$*-¥$2 + A, 

-269 =$8 — 8^7 -t-2lÇ5_20$3-f- MÇ^ 

26,2 = 512 — ii$io-+-45$8—84$6-i- 70$* — 2152+ H. 

M. Harzer a calculé des Tables qui donnent les logarithmes 
des 6 (jusqu'à 0<2) pour tous les centièmes de l'argument x, 
depuis x = o jusqu'à x = i {Astj\ Nachr,, 2759). Il a aussi 
donné deux exemples complètement développés, où l'application 
des formules est facilitée par l'emploi d'un procédé graphique 
{y^oiT BiilL, t. III, p. 61 5). 
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C'est an même ordre d'idées que se rattachent les recherches 
de M. Tchéb^chef sur les formules de quadrature à coefficients 
égaux, communiquées à l'Association française pour l'avancement 
des Sciences, à l'occasion du Congrès de Lyon [elles ont été 
reproduites dans le Journal de Mathématiques, 1878; voir aussi 
y Elude sur les formules d' approximation {ihid,, 1880)]. 

IV. 

45. Nous terminerons ces études par quelques remarques sur 
les formules de quadrature, qui sont étroitement liées aux for- 
mules d'interpolation." Il faut, pour cela, reprendre la notation 
symbolique du n" 8, en la complétant et la modifiant un peu. 

Désignant toujours par A les différences ordinaires, nous savons 
que l'opération (i -f- A)-'* a pour effet de changer j^o ^^y±ji' Quel- 
ques auteurs (Boole, Merrifield) remplacent i -|- A par le sym- 
bole E, de sorte qu'on a 

E = I -f- A, ^1 = E70, y±^n = E=^«J^o. 

Si nous réservons, d'autre part, le signe D pour les différences 
médianes, nous pourrons^ pour dédoubler l'opération D-, sup- 
poser qu'entre les termes donnés on ait intercalé ceux qui corres- 
pondent aux arguments fractionnaires db ^, di |, . . • et poser 

D2= =E + E-i — 2, 



(1) \ D = -— = =E2_E '2, 2v/i-HiD2 = E2-t-E «, 



L'opération fictive D se trouve ainsi parfaitement définie; mais, 
en fait, elle est remplacée par une moyenne arithmétique que, 
pour plus de clarté, je désignerai maintenant par la lettre M (et 
non plus par D, comme au n" 9), de sorte qu'on aura 

2M rr^AH ^ =E — E-1, 

1 -h A 



(2) { ^I =Ds/n-iD2, 

A = { D2 -f- M, 
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Le facteur y/i H-ID*-^ sert donc à passer de l'opération fictive D 
à l'opération réelle M, par laquelle s'obtiennent les différences 
intermédiaires; et si l'on veut confondre les différences d'ordre 
pair ou impair dans la même notation, et mettre, par exemple^ 
j)2/i+< au lieu de MD^'*, il n'y aura qu'à supprimer dans les déve- 
loppements le facteur y/i -\- jD-, ou à diviser par y/i +|D^. 

Il faut seulement distinguer entre les différences D^(o), relatives 

àj^oî et les différences D^(^) qui sont rapportées àyt_^= y/i + A.j^o* 
Les premières résultent des opérations 

I, M, D2, MD2, D*, ..., 
et les secondes des opérations 

^, D, MD, D3, MD3, .... 

Pour ramener les unes et les autres à la notation D^, D* , D^^ D^, 

D% . . . , il faut, dans le cas de l'argument (o), diviser par ^i + 1 D- 
les développements où ne figurent que les différences d'ordre impair 
et, dans le cas de l'argument (^), ceux où figurent les différences 
d'ordre pair. On a d'ailleurs 

r=i- :^D2-i- -i- D^ —, D6-f..... 



46. La différence A correspondant à l'accroissement fini h de 
l'argument x, on a, entre différences et dérivées, les relations 
symboliques 

( E = T -h A = e**, 

(3) 

^ ^ ( d=log(i + A), 

où, pour abréger, nous mettons d pour h -r-' On aura donc, en 
indiquant par des exposants négatifs des sommes ou des intégrales 

(4) 

( d^«=log^«(i4- A). 

Il est entendu que les exponentielles et les logarithmes doivent 
être remplacés parleurs développements en séries. Lagrange donne 
les coefficients de ces développements dans son Mémoire de 177'^ 
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Sur une nouvelle espèce de calcul {Œuvres, t. [II, p. 458). Dans 
les cas les plus simples, on a 

A =e*» — I =d-f-.î-d2-i- id3-f- -1 d*-4-..., 

2 6 24 

d.A-i= -^=i_id4--i-d5 Ld*-+.-^d6— ... 

e^ — I 1 12 7'20 30240 

— (-i)«7-^dî«H-..., 
^ ' {in)\ 

où B« est le n*^™® nombre de BernouUi ; pais 

d =log(TH-A) = A— - Aî+^ A3— i A^-h..., 

2 o 4 

Ad-i = î / ^, =h_Ia— — A2-4--Î-A3 

log(n-A) 2 12 24 

ï9 ., 3 ,, 863 ,^ 275 ,, 
720 160 60480 24192 

Ces développements permettent d'établir les formules usuelles de 
la manière la plus directe, en faisant usage des définitions 



d = A -^ , d-i 
dx 



- ,,J •••^^» 



n 



Le symbole A~% comme le remarque Boole, a un caractère inter- 
rogatif : la réponse consiste à renverser l'opération. Mais l'inté- 
grale finie A~* ou S n'est déterminée que par le point d'arrivée 
{yn^\)] le point de départ reste arbitraire; on est libre d'ajouter 

une constante d'intégration C, comme on fait pour l'intégrale /. 

L'indétermination cesse dès qu'on définit^ et / en indiquant la 

limite inférieure. 

Pour abréger les formules, je désignerai encore par S la somme 



n 




AS représente alors la somme des trapèzes contenus entre les 
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ordonnées j^OîJ^n • • «^JK/o q^^i figure dans la formule de quadra- 
ture la plus simple 

] y dx = hS. 



47. Les relations équivalentes 
(7) <i-^' r 



e^-i 



conduisent immédiatement à la formule sommatoire de Maclaurin 
ou d'Euler (') sous ses deux formes. La première donne, en dési- 
gnant par j-',^', ... les dérivées ^> ~Tî^ " ' ' 

ou bien 

ce qu'on peut écrire 

(8) ri-ro = -(yi"HX)--(rî-yo)+ — (ri^-n^)--... 

La seconde donne < 

ou bien 



iyVrf^=2^+ij.-Ay+ 



A3 

y — 



720 

et l'on voit que la formule d'Euler est à deux fins : elle peut servir 
de formule sommatoire (pour trouver 7^ j ou de formule de qua- 
drature (pour trouver j ydx\. Mais il faut la compléter, soit en 
ajoutant une constante d'intégration C, soit en introduisant deux 



(*) On peut consulter, sur cette formule, un Mémoire de Malmsten (1847), o" 
bien J. Tannery, Introduction à la Théorie des fonctions d'une variable, 
p. '352-363. 
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limiles, ce qui donne 

(9) lf"''ydx=S-^^ (^;_yj+^(y»_^™j— .., 

S étant la somme définie par l'expression (5). Le premier terme 
de la formule d'Euler est donc le second membre de la formule 
des trapèzes. 

Si, au lieu des relations (7), nous avions utilisé la suivante 

nous aurions trouvé la formule de Boole 

(10) A^ = J(y,+yo)-2(7'i+y;)+^(r^i+7î)-.-., 

puisque (2 + ^)yo=y^ -hyo- 

48. En partant de la relation 

(11) d-i=^ — -r^ — - =A-iH-- — — A -+-..., 

^ ' log(n-A) 2 12 

Lagrange obtient directement la formule de quadrature 

(12) l.J'ydx = j:y+ly-^\y^±^^y-.... 

Laplace (Méc. céL, t. IV, p. 206) y arrive par un détour. Il com- 
mence par démontrer la relation 

qu'il obtient en intégrant la formule d'interpolation de Newton, et 
qui s'écrirait directement 

log(i-hA) 

puis il en déduit la formule (12) par une sommation (SA = i). En 
introduisant des limites, cette formule devient 

(lO If r^^=S-~(A7«-A^o)+^(A2^«-Aî^o)-.... 
Mais les différences Ay,,, A-^«, . . . sont incommodes parce 
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qu'elles renferment des ordonnées situées en dehors des limites 
d'intégration; il y a lieu de recourir à une transformation par dif- 
férences ascendantes, en posant 

Iog(i-+-A) = — log^i- y^r^)' 
Il en résulte 

A [^ 2 I-+-A 12(H-A)2 24 (i-t-A)3 '"] 
— A-i I i A I A* 

~ 2 ~~ T2 T+~Â 24(l-+-A)2 '**' 

et l'on voit que les deux premiers termes sont les mêmes que 
dans (11), mais que les suivants sont transformés de façon que A^ 
devient =b A^_j. On trouve ainsi 

I r"^ I I 

C'est la formule que Laplace recommande comme « la plus 
simple et la plus commode ». 

49. Jusqu'ici nous avons fait usage des différences ordinaires A. 
Pour introduire les différences médianes D, on peut se servir des 
relations (i), (2), (3); mais, pour unifier la notation, il faut diviser 

par y/i 4- { D^ les développements où figurent D^^* (o) ou D-*(^). 
En ne considérant que. les opérations D, on aurait 

1 . __i . 
D = 62 — e i =2Sh|d, 



1 . 1 



(16) { 2 V^n- l D2 = e2 V e'^'^ = 2Gh I d, 

d = log(i -H I D2 + M) = 2 log ( J D + v/i-f-iD2), 

\d = argSh |D = argCh v/îTfD^. 

S'il s'agit des différences successives, il faut distinguer entre D" (o) 
et D''(^). On aura 

id*2« = [2argSh|D]^2«, 
v/n-iD2 

d=^2« =, ' [2 argSh J D]^^-", 

(i8)arg(4) { v/i + iD2 

di(2«+i) =: [2 argSh 1^ D]^(2«+i). 
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On retrouve ainsi, sous une autre forme, les relations très 
simples que M. Tisserand a indiquées dans une Note Sur un point 
du calcul des différences {Comptes rendus, 28 mars 1870). Les 
développements s'exécutent sans difficulté. Les formules connues, 
relatives aux fonctions hyperboliques, donnent 

2argShiD=D--^D3-i--LD5--4T:D7-H..., 
° ^ 24 640 7168 

(2argShiD)2 = D^ — — D*-h — D^— -^ Ds-h. . ., 

^ ^ ^ ^ 12 90 56o 

^ "o ^ ' 24 5760 967680 

2 argSh|_D ^P_.p3^ I p,_. p,_^ 
S/n-iD» 6 3o 140 

<l^I|M£ti = D-. - -L D + il D3 - J^ m -H . . . , 
4/, _4_iD2 '^ 7^^ 60480 

et ainsi de suite. On tire de là, par exemple, 

A^(o) = D(o)-|Da(o)+^D»(o)-..., 
ij 7t/^ = D-i(/i-hi)-D-i(i)+^[D(/i + i)-D(i)]-..., 



nh 



^J ^c/:p=D-i(7i)-D-i(o)-;l[D(/i)-D(o)]+.... 

Les D à exposants négatifs sont des sommes qui s'obtiennent en 
prolongeant le tableau des différences vers la gauche; ce sont les 
quantités désignées par '/} "ft • • • dans les Mémoires d'Encke sur 
les quadratures mécaniques (Jahrbuch, 1837 et 1862), où ce 
sujet est traité avec tous les détails qu'il comporte. 

50. On arrive aux mêmes résultats en partant des formules de 
Stirling (n'' 8), développées suivant les puissances de ^, la pre- 
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mière pour l'argument (o) 

-f- — ^ ( D3- i D5-H -^ D^-.. . ) 
et la seconde pour l'argument {{) 

» a \ 24 640 / 

(20) ( + - (d2- a D*+ #$■ D«-...)+ -^ fD3_ 1 D«+...) 

J 2 \ 24 5760 / 2.3\ 8 / 

_i!_ /d*- -^ D64-. . .W -^^ (D5-. ..) + .. .. 
2. j.4 \ 24 / 2.3.4.5 

Les coefficients des puissances de t représentent évidemment des 
dérivées. 

Nous écrirons maintenant /i à la place de t. 

Pour obtenir les expressions de d-', relatives k yn=/{nh)^ 
que nous désignerons par d-'(/i), il suffit de differentier ou 
d'intégrer les formules de Stirling par rapport à n. Nous pour- 
rions également les établir à l'aide des notations symboliques, en 
remarquant que 

/ d(/i) = E«d(o)=e«*ïd(o), 

d(/i)-4-d(— /i) A A/ \ 
^ = Gh/id.d(o), 



2 

(21) { d(/i)— d(— /i) 



= Sh/id.d(o), 



D = 2Sh-, ]\I = Shd. 
2 



Les fonctions hyperboliques Ch/zd et Sh/id étant exprimées par 



Shd et Sli -) on trouve d'abord 
2 



Chnô = i-\ D* H ^^ — — — ^ D* + 

1.2 1.2.3.4 



Shnd = M AH-Ai D^-hn ., , ^ D*-H... . 

L 1.2.3 1.2. 3. 4.5 J 
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On a ensuite, en divisant par le facteur de réduction y/i + {-D'^ la 
série des différences d'ordre impair, 

d(/i) + d(-/i) _ /^ ^ /i* _. \ 2ArgSh|D 



= H D* 



...) 



'^ \ 1.2 / y/i^lD2 

puis, en omettant le facteur de réduction dans le cas des diffc- 
rences d'ordre pair, et remplaçant M par ^ - ^ 

\ 6 12/ \i2o 36 90/ 

où les D se rapportent à l'argument (o). 

Nous n'insisterons pas davantage sur la démonstration des for- 
mules de ce genre, que l'on trouve longuement développées dans 
les Mémoires d'Encke et de Hansen {Relationen zwisclien Sum- 
men und Differenzen), Oppoizer, à son tour, a traité ce sujet 
dans le premier Chapitre du tome II de son Traité des Orbites 
{Différentiation et intégration numériques). Les formules pren- 
nent beaucoup de place, et elles sont d'ailleurs suffisamment con- 
nues. 

52. Des formules de quadrature d'un aspect très différent sont 
celles où figurent, non plus les différences successives des ordon- 
nées, mais les ordonnées elles-mêmes, et dont la plus simple est 
la formule des trapèzes. A cette catégorie appartiennent les for- 
mules de Cotes, qui supposent des ordonnées équidistantes ; 
celles de Gauss, qui sont moins commodes, mais qui permettent 
d'atteindre une précision double^ celles de Jacobi et celles de 
M. Tchébjchef, à coefficients égaux, etc. En dehors des Mémoires 
de Gauss, de Jacobi, AeT!c\iéhyc\\Qi {Journal de Mathé m,, 1873), 
on peut consulter ceux de Clausen et de Minding {Journal de 
Crelle, t. VI), Christoffel (/ètrf., t. LV), Mehler (/6tV/., t. LXIII), 
Encke {Jahrbuch, i863), Schellbach (1877), ^^ "^ travail plus 
récent de M. B. Baillaud (1881). J'ai résumé la plupart de ces tra- 
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vaux dans mon Etude sur les formules d^ approximation qui 
servent à calculer la valeur numérique d'une intégrale définie 
{Journal de Matliém., 1880; voir aussi Comptes rendus, 8 mars 
et 19 avril 1880, 16 juillet i883, et une Note de M. Callandreau, 
ibid.y 3 mai 1880). Il faut enfin citer les recherches de M. Stieltjes 
(Ann. de l'École Norm., 1884), de M. A. Markoff (/J/aM. Ann., 
1884), de M. J. Derujts {BulL de l'Acad. de Belg., 1886). On 
trouvera aussi quelques remarques intéressantes dans le Rapport 
de M. Merrifield. 

Je voudrais seulement, pour terminer, indiquer en quelques 
mots le parti qu'on peut tirer des formules (24) et (aS) du n^ 6, 
en supposant les ordonnées distribuées d'une manière symétrique. 
Soient donc y^^ j^^ y\, ^21 JK07 • • • les ordonnées correspondant 
aux arguments o, ±: a, li: 6, . . . , et posons 

yO = ^Oy = ^i, — <T2) • • • » 

1 1 

ri - y\ _> yt—y'i _^ 

T7. — ^1» TT — *^î» •••• 

Les formules en question pourront s'écrire 

y z= <Jq-\-Ox.X -\- <:'q,x^-\- 81 ,x{x^— a2) 



et 



^ = (Ti-t- Oi.a? -f- ci', {x^ — a2)-h Oj ,x{x^ — a-) 

^ aj {x^ — «2) (ar2 — 62 ) + o'[,x{x^ — a^){x^ — b^)^ . . , , 



ou 



•s "N (S <\ 

^, _ Qg — Qi _ Oi Qa 

^1- 62_a2 - «2- 62 "^ 62 — «2^ 



0« = 



•^z "^z > > ^ 

Oj — Oj Ôi Oj O3 



1 c2 — a2 (a2— 62)(a2_c2) (62 — a2)(62 — c2) ^c2 — a2)(c2 — 62) 
et d'une manière analogue 

, <Sy CTo , ^2 ^i 

^// __ ^J — ^O „„ _ ^ 2~^1 

et ainsi de suite. 
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L'intégration donne les formules de quadrature pour intervalles 



inégaux 



■^* . I . /i «2\ .. /i a2-f-62 «2^2 



i£ r^^ = -o+5^i-H(i-Ç).SH-(i- 






et 

dont la première suppose un nombre impair, et la seconde un 
nombre pair d'ordonnées. Pour avoir des intervalles égaux, il 
suffît de faire, dans la première formule, b = 2a, c = ia, . . ., 
et dans la seconde b = ia, c ^= 5a, .... En faisant encore, dans 

la première, a = — pour 2/2 + 1 ordonnées, et dans la seconde, 

a = pour 2 n ordonnées, de sorte que les ordonnées se trou- 
vent distribuées d'une manière uniforme entre les limites de l'inté- 
tégrale, on obtient les formules de Cotes. Si, enfin, nous prenons 
pour a, 6, c, ... les racines de l'équation P„ = o, où P„ est un poly- 
nôme de Legendre, nous avons la formule de quadrature de Gauss, 
car alors, en vertu d'un théorème connu, un certain nombre de 
coefficients de la formule s'annulent, et il en résulte que le degré 
de précision 2n — i s'obtient avec n ordonnées. 

53. Avec les coordonnées équidistantes j^o? y\ ? y\ > ^2» y^^ • • • 
correspondant aux arguments o, ± a, ± ia, ... (dont aucun ne 
coïncide avec 271), on peut encore construire une formule d'inter- 
polation qui mérite d'être mentionnée. Elle se compose d'expres- 
sions analogues aux fonctions symétriques D" du n** 2, où Xq, 
Xij X2, Xs, ... sont remplacés par les cosinus i, cosa, cos2ûr, 
cos3a, ... de sorte qu'elles prennent la forme 

DJ(/)= '^' 



(i — cosa)(i — cos:îa)(i — cos3a). . . 

H /l H.. 

{cosa — i)(cosa — cos*2a)(cosa — cos3a) . . . 

Les numérateurs /sont, d'une part, les quantités ^o, s^, S2y ... et 
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^i^ hf • • • qui dépendent des constantes données, 



2 

, _ ri — ri 

2 sina 



^2 = 



r^-^ri 



. _ r2— r2 

tj — : > 

2 sin2a 



. . . , 



et de l'autre les quantités i , cos^, cos2^, ... et — 
On a donc, par exemple, 



sin:r sin2â7 



sma sin2a 






^0 — ^1 
1 — cosa 



I — cosir 

Do(cosar)= , 



I — cosa 



D,(0 = 



ti t^ 



cosa — cos2a 



Si 



DJ(cosa7) = 



(ï — cosa)(i — cos2a) (cosa — i)(cosa — cos2a) 
^2 

(cos2a — i)(cos2a — cosa)' 

i cosa? 

_^ 

(I — cosa)(i — cos2a) (cosa — i)(cosa — cos2a) 

COS237 



DKO = 



(cos2a — i)(cos2a — cosa) 



+ 



h 



(cosa — cos2a)(cosa — cos3a) (cos2a — cosa)(cos2a — cos3a) 

h ^ 

(cos3a — cosa)(cos3a — cos^a) 

et ainsi de suite. On trouve aussi 



rDî?(cosar) 

(i — cosa7)(cosa — cosiPXcos2a — cosa?). . . [cos(n — i)a — cosrr] 
( I — cos na)(cos a — cos /ia)(cos 2 a — cos na) . . . [cos {n — i )a — cos naj 

2/*-i * \siQa/ 

__ sina7(cosa — cosa?)(cos2a — cosa?). . .[cos(/i — i)a — cosir] 
"~ sin na (cos a — cos /ia)(cos 2 a — cos na) . . . [cos {n — i ) a — cos na\ 

de sorte que DJ(cos;r) s'annule pour .r = o, a, 2a, , . . , (n — i)a 
et se réduit à 2"""* pour œ = na. Avec ces notations, le terme 
général de la formule d'interpolation peut s'écrire 



(-i)« 

2 



— -• (i — cosna){cosa — cosna). . .[cos(/i — i)a — cos/ia]DJ(5)DJ(cosa7) 
^~-sin/ia(i — cosna). . .[cos(/i — i)a -- cos/ia]Dî~*(/)Dj~* ( -; — j 
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et, en simplifiant l'expression des premiers termes, elle devient 



y = SQ-{-{Si—So) 



1 — cosa? 

■ ■ » * 

1 — cosa 
i(cos.a-,)(cosaa-cosa)D?(.)D5(cos..) 

- (cos 3 a — I ) (cos 3 a — cos a)(cos 3a — cos 2 a) DJ (5) D J (cos a?) -4- . . . 

1 . U — ti /sina? sinsa^X 
iSiniPH — sin2a ( -: ; ) 

2 cosa — cosaa \sina sin2a/ 

(SI II ayX 
-: 1-1-. . .. 
sina/ 

On a ainsi le développement de y suivant les sinus et cosinus 
des multiples de .x. 



FIN. 



ERRATA. 



Page 9, ligne 3 en remontant, au lieu de y^, lisez y^. 

Page 3i, note, ajoutez : Voir aussi les Remarques sur la théorie des moindres 
carrésj par E. Catalan {Mém. de l'Acad. de Belgique, 1878). 
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